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Abstract: We compute the one-loop scalar massless pentagon integral /f"^*^ in Z? = 6 — 2e 
dimensions in the limit of multi-Regge kinematics. This integral first contributes to the 
parity-odd part of the one- loop = 4 five-point MHV amplitude m'p at 0{e). In the 
high energy limit defined by s » si, S2 ^ —ti,—t2, the pentagon integral reduces to 
double sums or equivalently two- fold Mellin-Barnes integrals. By determining the 0{e) 
contribution to one therefore gains knowledge of m^^ through to 0{e'^) which is 

necessary for studies of the iterative structure of = 4 SYM amplitudes beyond one- 
loop. One immediate application is the extraction of the one-loop gluon-production vertex 
through to O(e^) and the iterative construction of the two-loop gluon-production vertex 
through to finite terms which is described in a companion paper [1]. The analytic methods 
we have used for evaluating the pentagon integral in the high energy limit may also be 
applied to the hexagon integral and may ultimately give information on the form of the 
R^^ remainder function. 
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1. Introduction 

In the planar M = 4 supersymmetric Yang-Mills theory, Anastasiou, Bern, Dixon and 
Kosower (ABDK) [2] proposed an iterative structure for the colour-stripped two-loop scat- 
tering amplitude with an arbitrary number n of external legs in a maximally-helicity violat- 
ing (MHV) configuration. The proposed iteration formula for the two-loop MHV amplitude 

(2)/ N • 

rUn [e) IS 

m'^He) = \ [m«(e)] ' + f^^\e) m«(2e) + C^^) + 0{e) , (1.1) 

thus the two-loop amplitude is determined in terms of the one-loop MHV amplitude rrin^ (e) 
evaluated through to 0{e^) in the dimensional-rcgularisation parameter e, a constant, C*-^-*, 
and a function, f^'^\e), which is related to the cusp [3, 4] and collinear [5, 6] anomalous 
dimensions. 

Subsequently, Bern, Dixon and Smirnov (BDS) proposed an all-loop resummation for- 
mula [7] for the colour-stripped n-point MHV amplitude, which implies a tower of iteration 
formulae, allowing one to determine the n-point amplitude at a given number of loops in 
terms of amplitudes with fewer loops, evaluated to higher orders of e. BDS checked that 
the ansatz was correct for the three-loop four-point amplitude, by evaluating analytically 
m^l\e) through to finite terms, as well as m^^\e) through to O(e^) and m^^\e) through 
to 0{e^). The BDS ansatz has been proven to be correct also for the two-loop five-point 

(2) 

amplitude [8, 9], for which (e) has been computed numerically through to finite terms, 
as well as 'm^c^\e) through to 0{e^). 

In the strong-coupling limit, Alday and Maldacena showed that the ansatz is violated 
for amplitudes with six or more legs [10]. This provoked the numerical calculation of 
mg^''(e) through to finite terms and of mg^''(e) through to 0{e^), where the BDS ansatz 
was demonstrated to fail in Ref. [11], and where it was shown that the finite pieces of the 
parity-even part of rn'^\e) are incorrectly determined by the ansatz^. In particular, it was 
shown that the remainder function, 

= (e) - \ (e)] ' - /(^) (e) m^^) (2e) - C^'^ , (1.2) 



^It was subsequently shown that the parity-odd part of mg (e) does satisfy the ansatz [12] 
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is different from zero for n = 6, where i?„ may be a function of the kincmatical parameters 
of the n-point amphtude, but a constant with respect to e. Because the calculation of 
Ref. [11] is numerical, the analytic form of Rq is unknown. 

Using the AdS/CFT correspondence, Alday and Maldacena [13] showed that in the 
strong-coupling limit planar scattering amplitudes exponentiate like in the BDS ansatz, and 
suggested that in the weak-coupling regime the vacuum expectation value of the n-edged 
Wilson loop could be related to the n-point MHV amplitude in = 4 super Yang-Mills 
theory. The agreement (up to a constant) between the light-like Wilson loop and the 
(parity-even part of the) MHV amplitude has been verified for the one-loop four-edged [14] 
and n-edged [15] Wilson loops, and for the two-loop four-edged [16], five-edged [17] and 
six-edged [18, 19] Wilson loops. Furthermore, it was shown that the light-like Wilson 
loop exhibits a conformal symmetry, and that the BDS ansatz is a solution of the ensuing 
Ward identities, up to functions of conformally invariant cross-ratios, which are present 
for n > 6 [17]. Recently, also the two-loop seven-edged and eight-edged Wilson loops 
have been computed numerically [20] , although no corresponding M = 4 MHV amplitudes 
are known. Thus, also the remainder functions r!j'^ and R^^ (for the Wilson loops, as 
well as for the MHV amplitudes, if we suppose that the agreement between them and the 
Wilson loops holds for seven or more points) are known numerically, and the numerical 
evidence [20] seems to confirm that they are functions of conformally invariant cross-ratios 
only. However, their analytic form is unknown. 

Through to finite terms, one-loop A/" = 4 MHV amplitudes, rrin^ (e), are parity-even and 
can be expressed in terms of box functions, which are massless for n = 4, one-mass forn = 5 
and two-mass easy box functions for n > 6 [21]. For n > 5, parity-odd contributions and 
one-loop n-edged polygons occur in the higher orders in e [22]. For example, an irreducible 
massless pentagon occurs in the parity-odd part of the one-loop five-point amplitude, which 
to all orders in e is [8] 

"^5'^ = -\^ .12523/^(1, 2, 3, 45, e) - |ei234/|"''(e) , (1.3) 

cyclic 

where m^^ denotes the one-loop coefficient, i.e., the one-loop amplitude rcscalcd by the 

tree-level amplitude, and where the cyclicity is over i = 1, ... ,5, /4™(1, 2, 3, 45, e) is the 

one-mass box with the massive leg of virtuality S45, J|~^^(e) is a pentagon evaluated in 

6 — 2e dimensions, and the parity-odd factor is £1234 = tr[75^i^2%/^4]- A massless hexagon 

occurs in the parity-even part [11, 22], as well as a pentagon in the parity-odd part [12, 22], 

of the one- loop six-point amplitude at 0{e). The hexagon may be reduced to one-mass 

pentagons. Thus, it seems fair to guess that the analytic form of the remainder function 
(2) 

Rq ' may be somewhat linked to the properties of the special functions occurring in the 
hexagon of the one-loop six-point amplitude. 

Because very little is known in the literature about the analytic properties of pentagons 

and hexagons to higher orders in e, in this paper we perform a first study of the pentagon 
occurring in the parity-odd part of the one-loop five-point amplitude (1.3). A massless 
pentagon can in general be reduced to quadruple sums dependent on the four independent 
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ratios of the kinematic invariants. In order to simplify matters, we choose a particular 
kinematic realm, the multi-Regge kinematics, where the quadruple sums reduce to double 
sums, which we proceed to study and expand to 0{e^). 

Our paper is organised as follows. As a prelude to our main objective which is the 
evaluation of the scalar massless pentagon integral in D = 6 — 2e dimensions in the multi- 
Regge limit, we first give a short review of multiple hypergeometric functions in Section 2. 
In general one may encounter a multiple hypergeometric function in several different rep- 
resentations - as a multiple sum, an Euler integral, a Mellin-Barnes representation or a 
Laplace integral. We then discuss in Section 3 different representations of Feynman inte- 
grals and show how they naturally match onto the various representations for hypergeo- 
metric functions. For example, in the Negative Dimension approach (NDIM) [23, 24] based 
on Schwinger parameters, one naturally finds the hypergeometric functions in the form of 
multiple sums. On the other hand, starting from Feynman parameters and utilising the 
Mellin-Barnes identity one naturally finds the hypergeometric functions as multiple Mellin- 
Barnes integrals. The one-loop pentagon integral and multi-Regge kinematics are defined 
in Section 4. In Euclidean kinematics, there are three distinct regions which we call I, 11(a) 
and 11(b). Regions 11(a) and 11(b) are related by symmetry. Regions I and 11(a) are related 
by analytic continuation. We use two distinct methods to compute the massless pentagon 
integral in D = 6 — 2e dimensions in each of these regions and show (a) that the solutions in 
each region are the same for the two approaches and (b) that the solutions in region I and 
11(a) are related by analytic continuation. In Section 5, we use NDIM to write the integral 
as multiple hypergeometric sums. In the multi-Regge kinematics these sums collapse to 
hypergeometric sums of two variables. It is straightforward to make a Laurent expansion 
around e = and find that hypergeometric sums can be written in terms of transcendental 
double sums which we denote by M functions. In Section 6, we revisit the problem start- 
ing from the Mellin-Barnes representation of the pentagon integral. In the multi-Regge 
kinematics we obtain two Mellin-Barnes integrals. We show that by summing the residues, 
one recovers exactly the same hypergeometric sums as in NDIM. Alternatively, by perform- 
ing the Mellin-Barnes integrals sequentially one can write the pentagon integral in terms 
of Goncharov's polylogarithms. In Section 7 we perform the analytic continuation of the 
pentagon to the physical region where all s-type invariants are positive. Appendices with 
more details concerning the generalised hypergeometric functions, nested harmonic sums, 
M functions and Goncharov's multiple polylogarithms we encounter in the calculation are 
enclosed. Further appendices contain some of the technical details omitted from Sections 5 
and 6 of the main text. 

2. Short review of hypergeometric functions 

In this section we briefly review the different representations of hypergeometric functions. 
We focus on Gauss' hypergeometric function, 

2F,ia,b,c;x) = Y^^-^^, (2.1) 

n=0 ^ 
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and for generalized hyper geometric functions we refer to Appendix C and to the litera- 
ture [25, 26, 27]. Eq. (2.1) defines the series representation of the hypergeometric function. 
The series is absolutely convergent inside the unit disc |a;| < 1. The hypergeometric func- 
tion however is defined over the complex plane, so we need to analytically continue the 
series (2.1) outside the unit disc according to the prescription x — > 1/x. Furthermore, in 
physical applications, one often encounters special cases where the indices a, b and c are of 
the form a + /3e, a, j3 being integers^, and it is desired to have the hypergeometric function 
in the form of a Laurent expansion up to a given order in e. We review in the following 
several other representations of the hypergeometric function, sometimes more suitable for 
specific problems. 

Let us first turn to the Euler integral representation, 

2F1 (a, 6, c; x) = ^) £dtt'-'{l- i)^"""^ (1 - xt)"". (2.2) 

This identity follows immediately from the fact that for < 1, we can insert the series 
expansion of (1 — xt)~"' in the integrand, 

y-(a)n n £(£) /"^^ b+n-w 

n! r{bnc-b)i 

^Y^Wn n r(c) r(5 + n)r(c-5) 

A. ril r(6)r(c-6) r(c + n) 

= 2-Fi(a,6, c;x), 

and the identity follows by analytic continuation for > 1. Note however that the Euler 
integral representation is meaningless if x is real and greater than 1 and a is a positive 
integer, since in that case the integral (2.2) is divergent. The Euler integral has the nice 
property that in the situation where all indices are of the form a+/3e, the expansion in e can 
be easily performed using integration-by-parts identities and the Laporta algorithm [28]. 
Indeed, since t and (1 — t) vanish at the integration limits, we can write, 

J^' dt ^ (t'-' (1 - ty-'-' (1 - xt)-'') = 0, (2.4) 

and carrying out the derivative on the integrand generates a set of recursive relations for 
the hypergeometric function. Using the Laporta algorithm we can solve the recursion and 
express every integral of Euler type in terms of a small set of master integrals. We can write 

down a set of differential equations for the master integrals that can be solved order by 
order in e in terms of harmonic poly logarithms (Sec Appendix B for a review of generalized 
polylogarithms) . Note that in the case of Gauss' hypergeometric function this procedure 

^In some applications, also half-integer values can be found. In the following we restrict the discussion 
to integer values. 
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(2.3) 



might look like an overkill, since we could as well expand the Pochhammer symbols in the 
series (2.1) and sum the resulting series in terms of nested sums. However, for more general 
hypergeometric functions, this naive approach might lead to series for which the sum is 
not necessarily known. 

A third way of representing a hypergeometric function is in terms of a Mellin-Barnes 
integral, 

^/ , X r(c) f+i°o r(a + z)r(b + z) , , , 

.F, (a, 6. X) = r(-.) r(c +l) (-^' ' ('■'^* 

where the contour is chosen according to the usual prescription, i.e. the contour should 
separate the poles in V{—z) from the poles in r(. . . + z). The identity (2.5) follows simply 
from the fact that for |x| < 1 we can close the integration contour to the left and sum up 
the residues of r(— z), which immediately results in the series (2.1). The Mellin-Barnes 
representation is however much more general, because the integral (2.5) is valid even for 
|x| > 1, in which case it suffices to close the contour to the right. In this sense, the 
integral (2.5) really represents the hypergeometric junction, because it is valid in all the 
regions^. This property gives us at the same time an effective way to perform the analytic 
continuation of the hypergeometric function outside the unit disc. Closing the integration 
contour to the right and taking residues at z = —a — n and z = — 6 — n, n G N, we find 



(2.6) 



2Fi(a, 6, c; x) = {-z)-" 2^1(0, 1 + a - c, 1 + a - 6; l/z) 

+ If^^^ ~ i-z)-" ,F,{b,l + b-c,l + b-a; l/z), 
1 (aji (c — b) 

where we used the formula 

r(a-n) = (-l)"-i^. (2.7) 
Finally, let us also mention the Laplace integral representation, 

T(c) 

2Fi(a,6,c;x) = -i-( / dti dt2 e-*^-*^ ij"^ ' (1 - *i *2 ^)-", (2.8) 
which follows directly from the integral representation of the F function, 

/•oo 

r{z)= / dte-H'-\ (2.9) 
Jo 

Since however in the following we do not use this representation explicitly, we do not 
comment on it further, but only quote it for completeness. 

The lesson to learn from this is that depending on what one wants to do, a given 
representation of the same hypergeometric function might be more convenient than another. 



^In Ref. [29], Appell and Kampe de Feriet propose for this reason to define the hypergeometric function 
through the MeUin- Barnes integral (2.5). 
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3. Feynman integrals 



In this section we review the different representations for n-point scalar one-loop Feynman 
integrals in D = Dq — 2e dimensions, Dq being a positive integer, 

= e-- / n (3-1) 

1=1 ' 

where the external momenta ki are incoming such that J2^=i ki = and the propagators 
have the form 



Di=k^- Ml + iO, 




(3.2) 



The external momentum scales are the Mandelstam variables Q|, and we work in the 
Euchdean region, Q? < 0. 

Feynman integrals can often be expressed in terms of (generalized) hypergeometric 
functions. In the previous section we showed that different representations of hypergeo- 
metric functions are useful to derive different properties. In the following we will argue that 
the different parametrizations used to evaluate Feynman integrals, Schwinger and Feyn- 
man parameters, series and Mellin-Barnes representations, are the equivalents to the four 
representations of the representations of the hypergeometric function, and switching from 
one parametrization to another might allow one to obtain valuable information about the 
Feynman integral. In the rest of this section we briefly review the different parametriza- 
tions, and in the rest of this work we give two examples how representations for Feynman 
integrals can be combined when computing Feynman integrals. 

3.1 The Schwinger parametrization 

The Schwinger parametrization^ is based on the identity 

Note that we explicitly derived Eq. (3.3) in Euclidean space, where we used the fact that 
in the Euclidean region Dj < in order to get a convergent integral. The corresponding 
relation in Minkowski space is similar, up to some factors of i in the exponent. Inserting 
Eq. (3.3) into the loop integral (3.1), we obtain, 

IS[{^^^}■,{Q^h{M^}) =/l?«/|^ exp(f^a,A), (3.4) 



where we introduced the shorthand 



Va = e^^^ n / da, a^'S (3.5) 



l\ m) Jo 



*Also known as a parametrization. 
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and performing the Gaussian integral in Eq. (3.4) leads to 

In{{^ih{Qlh{Mi}) = jva^ eMQ/V) eM-M). (3.6) 

The quantities "P, Q and ^A are polynomials in the Schwinger parameters a^, the internal 
masses Mf and the momentum scales Q?, 

n 

V = ^ai, 

i=l 

n-1 n /^~^ \ 

i=i j=i+i \ i=i / 

n 
i=l 

Note that the polynomials can directly be read off from the Feynman diagram in terms of 
trees and two-trees [30, 31, 32, 33]. 

The Schwinger parametrization is the starting point of the NDIM method, which we 
describe in the next section. Note however at this point the formal similarity between the 
Schwinger parametrization (3.7) and the Laplace integral representation of the hypergeo- 
metric function, Eq. (2.8). 

3.2 The Negative Dimension approach 

The crucial point in the NDIM approach is that the Gaussian integral (3.4) is an analytic 
function of the space-time dimension. Hence it is possible to consider D <Q and to make 
the definition [23, 24] 

"^""^ {er = n\5^^o. (3.8) 



for positive values of n. 

For the one-loop integrals we are interested in here, we view Eqs. (3.4) and (3.6) as 
existing in negative dimensions. Making the same series expansion of the exponential as 
above, Eq. (3.4) becomes 

dPk ^ {xiBiY^ 



ni,...,n„=0 1=1 



n rii 

Va /n(-ni,---,-«n;{Q?},{M2})[]^, (3.9) 

ni,...,n„=0 i=l 

where the are positive integers. The target loop integral is an infinite sum of (integrals 
over the Schwinger parameters of) loop integrals with negative powers of the propagators. 
Likewise, we expand the exponentials in Eq. (3.6) 



D OO 



(M;{Qf};{M.}) = / poE^^^^ E (^-lo) 



n=0 m=0 
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and introduce the integers qi, . . . ,qq, pi, . . . ,pn and mi, . . . , m„ to make multinomial ex- 
pansions of Q, V and M. respectively 



oo 


Qf 
gi! 




{qi + -- 




OO 

pi,...,p„=0 


Pi! 




(Pi + • • 




oo 

{-Mr = E 

mi,...,m„= 


(- 




mi ^_ 




=0 


mil 







(3.11) 

(mi + ... + mn)!, 



subject to the constraints 





= n, 


i=l 




n 






= —n — 


i=l 




n 






= m, 


i=l 






D 


..+qq 


~ ~ 2"' 



D 

2"' 



(3.12) 



Altogether, Eqs. (3.10) and (3.11) give 

I^{{M■dQn■,{M^}) = 



I 



Ua > j j j j j . . . j [pi + ...+p. 

Pn— U ^ 



Pi, 

qi,...,qq=0 
mi,...,mn=0 



(3.13) 



with the constraints expressed by Eq. (3.12). 

We recall that each of the Qi is a bilinear in the Schwinger parameters, so that the 
target loop integral is now an infinite sum of powers of the scales of the process (with 
each of the and the raised to a different summation variable) integrated over the 
Schwinger parameters. 

Equations (3.9) and (3.13) are two different expressions for the same quantity: 
Matching up powers of the Schwinger parameters, we obtain an expression for the loop 
integral with negative powers of the propagators in negative dimensions 

oo 

Pl,---,Pn—0 
qi,...,qq^O 



iir(i + m,)r(i + p.)J l^i_ir(i + ,,)i \ ^ 



n 

Pk 

k=l 



(3.14) 
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subject to the constraints Eq. (3.12). This is the main result of the negative dimension 
approach. The loop integral is written directly as an infinite sum. Given that Q can be read 
off directly from the Feynman graph, so can the precise form of Eq. (3.14) as well as the 
system of constraints. Of course, strictly speaking we have assumed that both Ui and D/2 
are negative integers and we must be careful in interpreting this result in the physically 
interesting domain where the i^i and D are all positive. Furthermore, of the many possible 
solutions, only those that converge in the appropriate kinematic region should be retained. 

3.3 The Feynman parametrization 

The Feynman parametrization is based on the identity 

n D-f^ - r{u,) . . . r(^„) l\ 7o {x,D, + ... + Xr^D^Y ' ^'-^'^ 

with V = Vi- Inserting this relation into Eq. (3.1), we obtain the following represen- 

tation for the Feynman integral, 

({^J; {Qh; (M.}) = rw j v. j (k^.^lk + jy ("^^ 

where we introduced the shorthand 
The coefficients Q and J are given by 

n i—1 



1=2 j=l 
n I i—1 \ n 

i=i \j=i / i=i 



(3.18) 



Performing the loop integration, we can reduce the Feynman integral to 

In ({^i}; {Qf }; {Mi}) = {-ir r{u - d/2) J Vx j^^, (3.19) 

where the F-polynomial is defined by 

F = Q^-J. (3.20) 



We have reduced the Feynman integral to an integral of a rational function over the unit 
cube in n dimensions. Note the formal equivalence of Eq. (3.19) with the Euler integral 
representation of the hyper geometric function, Eq. (2.2). 
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3.4 The Mellin-Barnes representation 

The MeUin-Barnes techniques rely on the following identity, 



1 



1 1 /•+*°° 



(A + B)* r(A) 

The contour in Eq. (3.21) is chosen in the standard way, i.e. it should separate the poles 
in T{—z) from the poles in r(A + z). We can apply Eq. (3.21) to the F-polynomial in 
Eq. (3.19), and break it up into monomials in the Feynman parameters Xi. The integration 
over the Feynman parameters can now be easily performed in terms of F functions, 



" r(ai)...r(a„) 



/ r[dxia:"'~^5(l -xi... -x„) = 
•^0 i=l 



F(ai + . . . + o„) 



(3.22) 



In this way we have eliminated all the Feynman parameter integrals in terms of Mellin- 
Barnes integrals, and we obtain a representation equivalent to the Mellin-Barnes represen- 
tation of the hypergeometric function, Eq. (2.5). 

4. The massless scalar pentagon in multi-Regge kinematics 

In the rest of this work we compute the scalar massless pentagon in D = 6 — 2e dimensions 
in multi-Regge kinematics. The scalar pentagon corresponds to the integral 

l5^(-i,-2,z.3,-4,-5;Q?) =e^- I ^D^D^BfD^^ ^^-^^ 

where the external momenta ki are lightlike, kf = 0, and are incoming so that Yl^=i ~ ^■ 
The external momentum scales are the Mandelstam variables = si2, S23, S45; si5 
with Sij = {ki + kj)'^, and we work in the Euclidean region, Sij < 0. Let us introduce the 
shorthands 



S = S12, Si =S45, S2 = S34, 
^1 = S51, t2 = S23- 



(4.2) 



We are concerned with the one-loop pentagon in the very peculiar kinematics, 

-S » -Si, -S2 > -ti, -t2- (4.3) 
The hierarchy of scales is such that 

S1S2 ~ sti ~ st2- (4.4) 
Equivalently, we can define the limit by the scaling 

S ^ S, Si-^Asi, S2^AS2, ti^X^ti, t2^>?t2, (4.5) 
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0.0 0.5 1.0 1.5 2.0 2.5 3.1 



Figure 1: The three regions contributing to the scalar massless pentagon in Euchdean kinematics. 



where A — > . The Euclidean region is itself divided into three other regions, in which 
the pentagon is represented by different analytic expressions. For later convenience let us 
introduce the following definitions, 



Sti ti St2 

X\ = = — and X2 = 

S1S2 K S1S2 

where we introduced the transverse momentum scale 



K ' 



— K = 



i-s) 



(4.6) 



(4.7) 



Note that k behaves as a t-type invariant under the scaling Eq. (4.5). In terms of these 
quantities the Euclidean region can be divided into three regions 

1. Region I, where ^/Xl + ^/x2 < 1. 

2. Region 11(a), where —^/xi + ^Jx^, > 1. 

3. Region 11(b), where — y^xi > 1. 

A graphical representation of these three regions in the {xi,X2) plane can be found in 
Fig. 1. Note that Region I is symmetric in xi and X2, whereas Regions 11(a) and 11(b) 
exchange their roles under an exchange of xi and X2- It is easy to see that Regions n(a) 
and 11(b) can be furthermore characterized by 

1. Region 11(a): {-h) < {-t2). 

2. Region 11(b): (-ti) > (-ts). 

Note that the region where ^4 is soft, si,S2 — ^ 0, corresponds to xi,X2 — ^ +00 in the 
(xi, a;2)-plane. 
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5. The pentagon integral from NDIM 
5.1 General considerations 

Solving the system of constraints from negative dimensions, we identify 125 quadruple 
scries contributing to the massless scalar pentagon in general kinematics. Each series has 
the form of a multiple generalized hypergeometric series. For example, 

j{ni,n2,n3,n4} ^ (_g)J'4B-f (^_f^Y51-§ ^_g^y345 + § ^_g^y45-§ (_i^)-^'612 + f 

# ^^^e r(^i)r(^2)r(i^3)r(z^4)r(z/5) 



-I T e 




r(z^345 - f )r(i^45i - f )rKi2 - #)r(f - i^45)r(f - u^i) 

D , D ^ D ^ D \ 

^45, y ~ ^51' ^ + ~ 2^345, 1 + y " ^451, ^ + ~ ^512; Xi,X2,X3,X4 I . 

(5.1) 

The arguments of the hypergeometric functions are ratios of scales, e.g. 

S2 SiS2 Siti ti 

Xl = —, X2 = — , X3 = ——, X4 = —, (5.2) 

S St2 St2 12 

and we introduced the definitions 1^123 = z^i + f 2 + ^3, etc. For convenience we have intro- 
duced the shorthand for quadruple sums, 

F{a, b, c,d, e, /; xi,X2, 2:3, ,14) 

(<2')ni+n2+n3+n4 (^)ni+n2+n3 {c)n2+n3+TH Xi^ X^ 2^3^ X^ (^-3) 



E 



ni,n2,n3,n4=0 (c^)m+n2 (e)n2+n3 (/)n3+n4 "l' "2! ^13! 714! 

The hierarchy of scales in multi-Regge kinematics eliminates many of the 125 solutions 
for the pentagon integral. The procedure for reducing the number of solutions is as follows, 

1. Any solution containing a summation that contains ratios of a "large" scale divided 
by a "small" scale, such as 

)n 
(5.4) 

cannot converge and is therefore discarded. This reduces the number of solutions 
from 125 to 22. 

2. Solutions with a pref actor that are less singular than 

111 , , 

SlS2 Stl St2 

when D = 6 — 2e and f i = 1 are discarded. This reduces the number of solutions 
from 22 to 20. 

3. Any sum that contains ratios of a "small" scale divided by a "large" scale such as 

-) , (5.6) 
s J 

gives its leading contribution when the summation variable n is zero. This leads to 
sums with fewer than four summations. 
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4. The remaining solutions contain only double sums of ratios of the three scales of 
Eq. (5.1), defined in Eq. (4.6). 

5.2 The solution in Region II 

As an example, we list here the solutions of the system of constraints in Region 11(a), 



St2 



with 



and 



^ND''^(s>Sl>S2,il,i2) = XI l'i^''\s,Si,S2M,t2), 



{I la). 



(5.7) 



(5.8) 



/{^^"^(s, si, S2, hM) = --^VT r(l - 2e) r(l + ef F4 (l - 2e, 1 - e, 1 - e, 1 - e; -yi, 2/2) , 
if ^"^(s, si,S2,hM) = ^ r(l + e) r(l - e) F4 (1, 1 - e, 1 - e, 1 + e; -2/1,2/2) , 



l'i^"'\s,Si,S2,ti,t2) = - 



r(l - 2e) r(-(5) r(-e - (5 + 1) r(e + (5) 



er(l + e)r(l-e)r(l-2e-(5) 



il-51-e-5 



l+el-e-5 1-5 



2/1,2/2 



4^^"^(s,Si,S2,tl,i2) = - 



r(-(5) r(l - 2e) r(-e - (5) r(e + 5) 



r(l + e) r(l - e) r(l - 2e - (5) 
X (-si)-^ 2/1+"^ 2/2"'"^ F4 (1, 1 - e, 1 + e + (5, 1 - e - ,5; -2/1 , 2/2) , 



(s,Sl,S2,tl,t2) 



er(l - e) r(5 + 1) r(-2e - (5 + 1) 



x(-si)-SjFo'2 



4'^'^"^(«'«l'*2,il,t2) = - 



1 1 + e 



1 



- 2/1,2/2 



1+e l+e+5 1-5 
V{-5) T{5 + 1) r(-e - 6f T{e + 5 + 1) r(l - 2e) 



r(l-e)2r(l + e)r(-2e-5 + l) 
X {-si)-^ 2/1+^ (1 + 5, 1 + e + (5, 1 + e + (5, 1 + e + (5; -2/1, 2/2) . 



where we defined 



Ik a^i ii 

yi = — = :r and 2/2 = — = — , 
a;2 t2 3^2 ^2 



and we factored out the usual loop factor 

r(i + 6)r(i-6)2 



rr 



r(l-2e) 



(5.9) 
(5.10) 

(5.11) 



The regulator 5 is needed in order to prevent divergences in the prcfactors in the limit 
Vi = 1 and D = Q — 2e. The singularities in 5 are spurious and cancel out in the sum over 



-14- 



all six contributions, e.g. , 



jijia) ^ j[iia) ^-^yl y-^ I [Inyi + V^(l - e) - V(-e)] ^4(1, 1 - e, 1 + e, 1 - e; -2/1,2/2) 



Mia) Alia) 
'fi T -la 



d 2,1 / 1 + 5 1 + <5 - e 



l + 1-e l + e + (5 



yi,2/2 



I (5=0 



[lnyi + V;(l + e)-V(-e)] ^4(1, 1 + e, 1 + e, 1 + e; -yi, ya) 



2 1 



1+6 l+6+e 



1 - - - 
1+S 1+e 1+e+S - 



- yi,y2 



1 5=0 



(5.12) 



The final result for the massless scalar pentagon in multi-Regge kinematics to all orders in 
e in Region 11(a) is then simply given by the sum 

= - ^ 2/2"' r(l - 2e) r(l + e)2 F4 (1 - 2e, 1 - e, 1 - e, 1 - e; -yi, ^2) 
+ ^ r(l + e) r(l - e) F4 (1, 1 - e, 1 - e, 1 + e; -yi, ^2) 



■ y^' < [ Inyi + V(l - e) - V'(-e)] F4 (l, 1 - e, 1 + e, 1 - e; -yi, ^2) 



^ ^2,1 I 1 + S 1 + S- e 
dS 



^ p>x ^0,2 



1+S 1-e 1+e+S - 



yuy2 



1 5=0 



+ ^yl\[\nyi+tp{l + e)- iP{-e)] ^4(1, 1 + e, 1 + e, 1 + e; -yi, y2) 



(5.13) 



^ 5 2,1 1 + S 1 + S + e 
^dS^^'^i - 



1+S 1+e 1+e+S - 



- yi,y2 



1 5=0 



Note that the only functional dependence of l^^i^ is in the ratio of scales y\ and y2) ^-e. , 
in the transverse momentum scales ii, ^2 and k, 



(5.14) 



The solution in Region 11(b) is related to the Region 11(a) by analytic continuation accord- 
ing to the prescription ti/t2 ^ ^2/^1, or equivalently y2 l/y2- From the symmetry of 
the multi-Regge limit in t\ and t2 it is easy to see that we must have 



11 



(5.15) 



In Appendix E we explicitly show that Eq. (5.13) enjoys this property. 
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Using the reduction formulas for the Appell function given in Appendix C, we could 
reexpress all the F4 in Eq. (5.13) in terms of Gauss' hypergeometric function, 

F4(l-26,l-6,l-e,l-e, n _ _ n _ _ J 



^4(1,1- 



(l-x)(l-y)'(l-x)(l-y). 
= (1 - x)^-' (1 - y)^-' 2^1(1 - 2e, 1 - e, 1 - e; xy) 
_ {l-x)^-'{l-y)^-' 
(1 - xy)^-^' 

e,l + e, 1 — e, 



(5.16) 



:i-z)(l~y)' (l-x){l-y), 
= (1 - x) (1 - y) Fi(l, 2e, 1 - e, 1 + e; xy). 

The Appell Fi function appearing in this reduction can be easily expanded into a Laurent 
series in e using XSummer [34]. Note however that we do not know the corresponding 
reduction formulas for the Kampe de Feriet function appearing in Eq. (5.13). For this 
reason, we do not apply the reduction formulas of the Appell F4 functions, but we proceed 
and perform the e expansion directly on the series representation of the hypergeometric 
functions appearing in Eq. (5.13). Since all the hypergeometric functions in Eq. (5.13) are 
finite for e = 0, wc can safely expand the Pochhammer symbols into a power series under 
the summation sign, 

(l + e)„ =n! {l + eZi{n) + e^Zuin) + e^Zui{n) + 0{e^)), 



= - (1 - e 5i (n) - (n) - (n) + 0{e^)) , 



(l + e)„ n! 

where S and Z denote nested harmonic sums and Euler-Zagier sums, defined recursively 
by [35] 

n 

S,{n) = Z,{n)=H(i^=^-, 

''^^ (5.18) 

k=l k=l 
Using the algorithm described in Appendix A, we can express all the Euler-Zagier sums 
in terms of harmonic sums. Finally, using the algebra properties of the S-sums, we can 

reduce all the products of harmonic sums to linear combinations of the latter. Inserting 
Eq. (5.17) into the series representation for hypergeometric functions we obtain the desired 
e expansions, e.g. , 

Fi{l, 1 + e, 1 + e, 1 + e; xi, X2) = M{0, 0, 0; .xi, ,12) 

+ e M{0,0,l;xi,X2) - A1(l, 0, 0; xi, X2) - A1(0, 1, 0; xi, X2) 

+ [X((l, 1), 0, 0; XI , X2) + M(0, (1, 1), 0; xi , X2) + MiO, 0, (1, 1); xi , X2) 

+ M{l,l,0;xi,X2) - A1(l,0, l;xi,X2) - A^(0, 1, 1; xi, X2) - >J(0, 0, 2; xi, X2) 
+ 0{e^). 

(5.19) 
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The Ai functions appearing in this expansion are transcendental functions defined by the 
double series 

OO OO ^ X 2 

Mii,j,k;x,,X2)= J2 E [n) Sririi) Sjin^) S^{m + x^' x^\ (5.20) 

ni=0n2=0 V / 

Note that due to the appearance of the binomial squared term in Eq. (5.20), we cannot 
reduce the double sums in general to known function using the standard techniques [34, 36]. 
We can however sum the series in some particular cases in which we can relate the M- 
function to the expansion of a known hypergeometric function. This issue will be addressed 
in Appendix D. 

Since Eq. (5.13) only involves AppcU functions and Kampc de Ferict functions with 
indices 1 + Qe, it can be easily expanded in terms of M. functions. The first two orders 
read, 

li;i''\K,h,t2) = ii'"'\yi,y2) + ei['"'\yi,y2) + Oie% (5.21) 

with 

4''"^ (yi, 2/2) = (-8 In yi - 4 In y2)M (O, 0, (1, 1); -yi , ya) - 4 In ysA^ ((1, 1), 0, 0; -yi , ys) 
+ 18M (0, 0, (1, 2); -yi , 2/2) + 18M (O, 0, (2, 1); -yi , 1/2) - 24M (0, 0, (1, 1, 1); -yi , ^2) 
+ 8M{0, 1, (1, 1); -yuy2) + 16Al(l,0, (1, 1); -yuy2) - 8M{{1, 1),0, 1; -2/1,2/2) 

+ ((1, 1), 1, 0; , ,2) - M (0, 0, 0; -y, , ^2) + '-^^ + ^ - 2(3) 

2 In yi In y2 + In^ y^ + — J + (6 In yi + 3 In y2)M (0, 0, 2; -yi , 1/2) 

+ (2 In 2/1 In y2 + ^) X (1, 0, 0; -yi , ^2) + (4 In yi + 4 In y2)M (l, 0, 1; -yi , ^2) 

+ 4 In yiM (0, 1, 1; -yi , y2) - 4 In yiM (l, 1, 0; -yi , y2) + ( In^ yi + tt^) M{0,1, 0; -yi , y2) 

+ In y2M (2, 0, 0; -yi , y2) - 12M (O, 0, 3; -yi, ys) - QM (0, 1, 2; -yi, ys) 

- 127W (1, 0, 2; -yi, ys) - 87W (l, 1, 1; -yi, ys) + 2M (2, 0, 1; -yi , ys) 
-27W(2,l,0;-yi,y2), 

(5.22) 

4''"^ iyuy2) = (0, 0, (1, 1) ; -yi , y2) (4 In yi In y2 - 4 In^ yi + 2 In^ y2 + Att^^ 

+ (2 ln2 y2 - 4 In yi In y2)A^ ((1, 1) , 0, 0; -yi , ys) + (8 In yi - 12 In y2)A^ (l, 0, (1, 1); -yi , y2) 
+ (4 In y2 - 8 In yi)Al ((1, 1), 0, 1; -yi , y2) + (8 In yi - 4 In y2)A^ ((1, 1), 1, 0; -yi , ys) 

- 15 In y2M (O, 0, (1, 2); -yi , 2/2) - 15 In y2A^ (O, 0, (2, 1); -yi , y2) 
+ 20 In y2A^ (O, 0, (1, 1, 1) : -yi , y2) - 4 In y2Al (0, 1, (1, 1); -yi , y2) 

- In y2A^ ((1, 2), 0, 0; -yi , y2) - In y27W ((2, 1), 0, 0; -yi , y2) 

(5.23) 
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+ 4 In ya-M ((1, 1, 1), 0, 0; -yi , yz) + 327W (0, 0, (1, 3) ; -yi , ya) + 36M (0, 0, (2, 2); -yi , ya) 
+ 327W (0, 0, (3, 1); -yi , ys) - 48A^ (O, 0, (1, 1, 2); -yi , ys) - ASM (0, 0, (1, 2, 1); -yi , ys) 

- 48M (0, 0, (2, 1, 1); -yi , ys) + (O, 0, (1, 1,1,1); -yi , ys) + 12A^ (0, 1, (1, 2); -yi , ya) 
+ 12M (0, 1, (2, 1); -yi, y2) - WM (0, 1, (1, 1, 1); -yi , ys) + 18A1 (l, 0, (1, 2); -yi, ys) 

+ 18M(1,0, (2, 1); -yi,y2) - 247W(l,0, (1, 1, 1); -yi,y2) + 8M{l, 1, (1, 1); -yi,y2) 

- 27W((1, 2), 0, 1; -yi, y2) + 2M ((1, 2), 1, 0; -yi, y2) - 2A1 ((2, 1), 0, 1; -yi , y2) 

+ 2>l((2,l),l,0;-yi,y2) +8M((l,l,l),0,l;-yi,y2) -8A^((l,l,l),l,0;-yi,y2) 

TT^ In yi In^ yi In y2 2 In^ yi Stt^ In y2 



+ (0,0,1; 


-yi 


y2) 


+ A1(0,1,0; 


-yi 


y2) 


+ X(1,0,0; 


-yi 


2/2) 


+ M(0,0,2; 


-yi 


y^) 


+ M(0,1,1; 


-yi 


2/2) 


+ 7W(1,1,0; 


-yi 


^2) 


+ A^(2,0,0; 


-yi 


^2) 



lnyilny2- 3 ^ 3 

TT^lnyi ln^yilny2 21n^yi 7r^lny2 



+ 



6C3) 



+ 



- Inyi In y2 + 



+ 2C3 

TT^lnyi 31n^yilny2 7r^lny2 



+ 



+ 2C3) 



3 In yi In y2 + 3 In^ yi 



3 In y2 



-37r2) 



-21nyi lny2 + 21n^ yi 

2 In yi In y2 - 3 In^ y^ + — j + (In y2 - 2 In yi)7W (2, 1, 0; -yi , y2) 



In yi In y2 



ln^y2 



+ (9 In y2 - 6 In yi)7W (1, 0, 2; -yi , y2) 



01 ^ , 7rMn^y2 , 2t,\ 
21ny2C3 + 1 + -^) 



+ (4 In y2 - 4 In yi)7W (1, 1, 1; -yi , y2) + (2 In yi - In y2)7W (2, 0, 1; -yi , y2) 

_L A/fCn n n x/ln2yiln2y2 tt^ In^ yi In3yilny2 1^2 
+ M (0, 0, 0; -yi , y2) 1 g 3 

+ (3 In^ yi - 2 In^ y2 - tt") (1, 0, 1; -yi , y2) + 10 In y2M (0, 0, 3; -yi , y2) 

+ 3 In y2A^ (0, 1, 2; -yi , y2) - (O, 0, 4; -yi , y2) - 87^ (0, 1, 3; -yi , y2) 
- 12^ (1, 0, 3; -yi , y2) - (l, 1, 2; -yi, y2) . 

Note that i^^"^ and 4^^"^ are of uniform transcendental weight 3 and 4, as expected. 

5.3 The solution in Region I 

The solutions in Region I read 



/^-^)(s,si,S2,ii,i2) = rve"'''^^^^ — l^l^{s,si,S2,tiM), 



S1S2 



(5.24) 



with 



2ND(^'«l'«2,il,t2) = ^(s,Sl,S2,il,i2), 



(5.25) 



1=1 
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and 



Il\s,Si,S2,ti,t2) 



x^'x2'T{l-2e)T{l + e) 



X ^4(1 — 2e, 1 — e, 1 — e, 1 — e; —xi, —X2), 
^ {s, si, S2MM = A r(l + e) r(l - e) ^4(1, 1 + e, 1 + e, 1 + e; -XI, -X2), 



eT{l + e) 



00 00 



Ml-".2 



ni=0 n2=0 



X 



(-e - (5 + l)ni ni!n2! 

* {-S2f' i-h)-' XY 

r(l - 25) T{5)T{-e -6 + l)T{5- e) T{e + 5) 



a;2 



r(i-5)r(,5 + i)r(i + e)r(i-e)2 



X F\ 



,3,0 



1,2 



1 1 - e - (5 1 - 25 

1-5 



l-€-5 l+e-5 



Xi, -X2 



l5\s,Si,S2,ti,t2) = -{-sf {-Si) (-^2)" X. 

T{-5) r{25 + 1) r(-e - 5) T{e - 5) T{-e + 5 + 1) 



X 



r(i-5)r(5 + i)r(i + e)r(i-e)2 

1 l-€ + 51 + 25 
- - l + 



l + e + (51-e + (5 



-Xi,-X2 



Ie\s,si,S2,h,t2) = - {-S2)^ (-si) 



r(5)Ti-e-5) 
er(l-e) 



00 00 



xEE 

ni=0 n2=0 



ni— 712 ((5 + 1)^2(1 -e)n2(e) 
(e + <5 + l)„,ni!n2! \x2j ^ ^' 



(5.26) 



The regulator 5 is introduced to prevent divergences in the T functions in the prefactor. 
The cancellation of the spurious 5-poles is not as straightforward as in Region 11(a), due the 
appearance of a new type of hypergeometric series. Since however this new series involves 
only Pochhammer symbols of the form (.)„^_„2, {■)ni and (.)n2) it can be reduced to Kampe 
de Feriet functions using the techniques described in Appendix C. We find 



00 00 

EE 

ni=0 712=0 



(1 - 5)n2 (-5)ni-n2 (<^ + l)n2 (1 " e)n2 ("C) 



Til 



(-^2) 



712 



xi 0,3 
l-e-5 X2 



{-€-6 + l)ni ni!n2! 

1 + 5 1 1-5 1-5 1-el-e 



2 2-€-5 



- Xl, 



X2 



1-5 1+5 1-e 



1 - - 
1+5 1-5-e 1+e 



Xl, -X2 



(5.27) 
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EE 

ni=0 n2=0 



n\—n2 



)n2 V^)n\—n2 



{e + 5 + l)nini\n2\ 



ni 



"2 



5( Xl_ 0,3 



2 2 + e + (5 



1 1 + 5 1 + 5 1-e 1 + e 



3^2 



l + 1-e 



l-(5 l + 5 + e 1-e 



(5.28) 



After inserting these expressions into Eq. (5.26), we expand the solution into a Laurent 
series in 5. The poles in 5 cancel mutually between /g^^ and I^^ and I^^ and Iq \ After 
some algebra, we find the following expression for the solution in Region I, 

= - ^ xj-" X2 ' r(l - 2e) r(l + ef Fi{l - 2e, 1 - e, 1 - e, 1 - e; -xi, -X2) 
+ ^ r(l - e)r(l + e) F4(l, 1 + e, 1 + e, 1 + e; -xi, -X2) 

— xj;"^ < [lnx2 + -0(1 — e) — V'(— e)] -F4(l, 1 — e, 1 — e, 1 + e; — xi, — X2) 



+ 



7?^ 

1-e X2 



2 2-e 



- - 1 
1-e l+e+5 -1+6 

1 1 1 1 1 - e 1 - e 



Xi,-X2 



1 5=0 



(5.29) 



- xi, ■ 



Xl 



2^2 



e2 



[lnx2 + ^'(l - e) - V'(e)] ^4(1, 1 - e, 1 + e, 1 - e; -xi, -X2) 



9 2,1 / 1 + (5 1 - e + 5 
+ a5^0'2 - 



e ^1 170,3 
— -'^2,0 



1 + e X2 



2 2 + e 



- - 1 
1+e l-e+5 -1+6 

1 1 1 1 1 - e 1 + e 



Xl, -X2 



|a=o 



Xl 
Xl, — 
X2 , 



Note that the right-hand side of Eq. (5.29) only depends on the dimensionless quantities 
Xl and X2, i.e. , on the transverse momentum scales k, ti and t2, 

I^l,{s,Si,S2,ti,t2) =I^^{K,ti,t2). (5.30) 

Furthermore, we know that the pentagon in Region I must fulfill the symmetry relation 

(5.31) 



2^Nd('^'*1'^2) =I^l{K,t2,ti). 



The solution given in Eq. (5.29) however apparently breaks this symmetry, due to the ap- 
pearance of the ratio Xi/x2- We show in Appendix E that Eq. (5.29) indeed has the correct 
symmetry properties, which becomes explicit only after a proper analytic continuation has 
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been performed. We also show that the solution in Region I, Eq. (5.29), can be obtained 
from the solution in Region n(a) by performing analytic continuation according to the 
prescription yi 

6. The pentagon integral from Mellin- Barnes integrals 

In Ref. [37] a Mellin-Barnes representation for the pentagon was given, 



T{D-N,) 

X r (z/5 + + 22) r - AT^ + Z/i - Zi - 22 - Z-^ r (^2 +Z2 + Z3) 

X r - Ni, + - Z2 - Z3 - r (1/4 + Z3 + Z4) 

xt(-^ + N^ + Zi+Z2 + Z3 + Z4], 



(6.1) 



where N^, = ^Vi- In particular, if we put Vi = 1 and I? = 6 — 2e, then we get 
75^(1,1, l,l,l;g2) = 



_p7Be (_s)-2-e 



r(l-2e) 

4 

I y ■ "™ . 

X 



X r (Zi + Z2 + 1) r (-e -ZI-Z2- Z3 - 1) r (Z2 + ^3 + 1) 
Xr{-e-Z2-Z3-Z4-l)T {Z3 + Z4 + l)r {e + Zi + Z2 + Z3 + Z4 + 2) . 

In the following we extract the leading behavior of this integral in the limit defined by the 
scaling (4.5), and we show that this leading behavior is described in all the regions by a 
twofold Mellin-Barnes integral, which can be evaluated in terms of multiple polylogarithms. 

6.1 Evaluation of the Mellin-Barnes integral in Region I 

Performing this rescaling (4.5) in the Mellin-Barnes representation (6.2) and after the 
change of variable Z4 = z — zi — Z2 — Z3, we find 



15^(1,1, l,l,l;Q,f) 
1 

X 7 r 

(27ri) 



r(l-2e) 

L,/_;;a.a...a.(^)~-(f)— (^)''(7)"- 

X T {-e - zi - l)T {-€ - z + zi - l)T {z - zi - Z2 + 1)T {-Z2) T {zi - Z3 + I) 

XT{€ + Z- Z2- Z3 + 2)T {-Z3) r {Z2 + Z3 + 1) 
X r {-Zi + Z2 + Z3) r {-Z + Zi+ Z2 + Z3) . 

(6.3) 
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To obtain the leading behavior in our hmit A ^ let us follow the strategy formulated, 
e.g., in Chap. 4 of [32, 33]. We think of the integration over z as the last one, and we 
analyze how poles in r(. . . — z) with leading behavior might arise. There is only 
one possibility, coming from the product F (— e — z + zi — 1) F (— e — z\ — \). Taking the 
residues at zi = — 1 — e + ni, ni G N, we find 

A^F (-e - z + zi - 1) ^ A^F (-2e - 2 - z + m) . (6.4) 

If we now take the residues at z = — 2 — 2e + ni + n2, n2 G N, we find 

A^F (-2e - 2 - Z + m) ^ _^-2-26+ni+n2 _ (-Q_5) 

Since we are only interested in the leading behavior in A~^, we only keep the terms in 
n\ = n2 = 0. Hence, we find a twofold Mellin-Barnes representation for the pentagon in 
multi-Regge kinematics, 

75^(1,1, 1, 1, 1; Qf ) = rr e^-^ iz^jW t^^f,) (6.6) 

with 

= r(i + .)"r(i-.)2 (2^*? U r (-.0 

X F (-e - Z2) F (-e - zi - Z2) F {-Z2) F {zi + Z2 + l)r {e + zi + Z2 + if , 

(6.7) 

where xi and X2 are defined in Eq. (4.6). Note that this expression is symmetric in xi 
and X2, as expected in Region I. We checked that if we close the integration contours to 
the right, and take residues, we reproduce exactly the expression of the pentagon obtained 
from NDIM, Eq. (5.29). 

We now evaluate the Mellin-Barnes representation (6.6) and we derive an Euler integral 
representation for the pentagon in multi-Regge kinematics. Let us concentrate only on the 
Mellin-Barnes integral. We start with the change of variable Z2 = z — zi and we find 

rn -1 1 r^'°° 

^-('^'^-^^) = rxTTW^(^y-.oo d.^d..r.r-r(-.-.)F(. + i) 

xT{e + z + if F (-e - zi) F {-zi) T {zi - z)T {-e - z + zi) . 

We now replace the Mellin-Barnes integral over zi by an Euler integral by using the trans- 
formation formula, 

f +ioo 

dzi r{-zi) r(c - zi) r(6 + zi) F(a + zi) X^i 



- / 



= F(a) F(6 + c) dv u^-i (1 - i;)«+^-i (1 - (1 - ^)t')~", 
^0 



(6.9) 



and we find 



llihK,ti,t2) =r- -^^ r:.— [^"^dz f\v(l-v(l-—\]^ X^ 

mbV , 1, 2; r(l + e) r(l - e)2 27ri y_,^ 70 V V ^2)) ' (6.10) 
X (1 - v)-^'-'-^ v-'-^ r(-e - zfT{z + l)r(e + z + if. 
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To continue, we exchange the Euler and the Mehin-Barnes integration. Note that we must 
be careful when doing this, because some of the poles in z could be generated by the Euler 
integration. We checked numerically that in the present case the exchange of the two 
integrations is allowed and produces the same answer. We now close the z-contour to the 
right and take residues at 2; = n — e, n G N. 



X { ln^(l-'y) + ln^v + ln^ ( \-v[ 1 " X2 ^ -^{n + \f ^ + n + \f 

+ 7r^ + 21n(l-i;)ln?;-21n(l-i;)ln^l-?; {\ 

- 21ni;ln {^-v{^ - — - 21n(l - t;)lnx2 - 21ni;lna;2 

+ 21n ^1 - u ^1 - j lnx2 + 21n(l - v) V'(n + 1) + 2\x^vi){n + 1) 

- 21n {\-v{\- 1) -21na;2V'(n + 1) - 21n(l - v) V'(-e + n + 1) 

- 21nt;V'(-e + n + 1) + 21n {\-v{\ - + " + 1) 
+ 2 In a;2 V'(-e + n + 1) - 2'0(n + l)'0(-e + n + 1) - (^^ + 1) 

(6.11) 

To continue, we perform the e-expansion under the integration sign^ 

X£^(«,ti,t2) =xi^)(xi,X2) + ejJ^)(xi,X2) +0(e2). (6.12) 



We find 



and 



iy{xi,X2) = dv^ ■ , (6.13) 

Jq — XlV + X2V — V — X2 



T\\xi,X2)= d^^^— 7 7-\ • ^^-^^^ 

Jo + [-Xi+ X2- l)v - X2 

where i^^^ and i^^^ are functions depending on (poly) logarithms of weight 2 and 3 respec- 
tively in xi , X2 and v (See Appendix F for the explicit expressions) . Note that this implies 



^Note that we must be careful when we do this, since the integrals might contain poles coming from the 
terms l/v{l — v). As we will see in the following however these poles are spurious and cancel out. 
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that Xq {xi,X2) and Xj (xi, ,^2) will have uniform weight 3 and 4 respectively, as expected. 
Furthermore note that the poles in v = and v = 1 have cancelled out. However, we still 
need to be careful with the quadratic polynomial in the denominator of the integrand, 
since it might vanish in the integration region. We analyze this situation in the rest of this 
section. 

We know already that the phase space boundaries in Region I require 

v/ir + v^<l. (6.15) 

This subspace of the square [0, 1] x [0, 1] is at the same time the domain of the integral 
X(xi, 2:2; e). We can divide this domain further into 

1. Region 1(a): xi < X2- 

2. Region 1(b): X2 < xi. 

We now turn to the quadratic denominator in Eqs. (6.13) and (6.14). The roots of this 
quadratic polynomial are 



Ai = Ai(xi, X2) = (1 + X1-X2- -s/A^) , 

I ; ! (6.16) 

A2 = Ai(a:i,X2) = - yl + xi-X2 + \/XkJ , 



where A^ denotes the Kallen function 

Xk = \k{xi,X2) = \{xi,X2, —1) = 1 + xl+ X2 + 2xi + 2X2 — 2X1X2- (6-17) 

First, let us note that Ax(a^i, 3^2) > 0, V (xi, X2) G [0, 1] x [0, 1], and hence the square root 
in Eq. (6.16) is well defined in the Region I. Second, it is easy to show that on the square 
[0, 1] X [0, 1] we have, 

-1 < Ai(xi,X2) < and 1 < A2(xi, X2) < 2. (6.18) 

For later convenience, let us note at this point the following useful identities 

A1A2 = -X2, 
Al + A2 =1+Xi — X2, 

Ai-A2=-V^, (6-19) 
Al) (""^ A2) X2' 

From Eq. (6.18) it follows now immediately that the quadratic denominators in Eqs. (6.13) 

and (6.14) do not vanish in the whole integration range [0, 1], and hence all the integrals 
in Eqs. (6.13) and (6.14) are convergent. Using partial fractioning and the relations (6.19) 
we can write 

X«(xi,X2) = ^ ^^ ^^'Hx„X2,v) _ 1 ^ J'^Hx„X2,v) ^ (g_20) 

V Ak Jo V- X2 VXk Jo v - Xi 
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and 



V Jo V- X2 vXk Jo v - Xi 



Let us introduce at this point the function 

A3 ^ Xsix,,x,) = = , ^ . (6.22) 

X2- Xi Ai + A2 - 1 

This function appears in i^^^ and i^^^ through the logarithm 

It is easy to see that this logarithm is well defined in Region I, since 

1. In Region 1(a), xi < X2, and hence A3 > 1. 

2. In Region 1(b), X2 < xi, and hence A3 < 0. 

Note however that A3 diverges on the diagonal xi = X2- This is not a contradiction since 

lini / -= hm / ^V=0, (6.24) 

X2-*xiJq t- X3{xi,X2) Xs-^ooJo t-Xs 

in agreement with 

lim In ( 1 - 1; ( 1 - — ) ) =0. (6.25) 



X2^xi Y \ X2 



We now evaluate the integrals Xq ' and Xj ' explicitly. To this effect, let us introduce 
some generalized multiple polylogarithms defined by 

G{a,w;z)= [ dt f{a,t) G{w;t), (6.26) 
Jo 



where 



If all indices are zero we define 



/(a,i) = ^. (6.27) 



dt ^ 1 
G(0„; z)= - G(0„_i; i) = - In" z. (6.28) 
Jit n. 



In particular cases the G-functions reduce to ordinary logarithms and polylogarithms, 

G(a„; z) = -^ In'* f 1 - -) , ^(On-i, a; z) = -Li„ (-) . (6.29) 



Note that these definitions are straightforward generalizations of the harmonic polyloga- 
rithms, and hence these functions inherit all the properties of the HPL's. In particular 
they fulfill a shuffle algebra 

G{wi;z)G{w2;z) = G{w;z). (6.30) 
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If the weight vector w has a parametric dependence on a second variable, then we obtain 
multidimensional harmonic polylogarithms. Finally, let us introduce the following set of 
functions, which will be useful to write down the answer for the pentagon 

M{w) = G{w;l). (6.31) 

The M-functions defined in this way are in fact Goncharov's multiple polylogarithm (up to a 
sign). For a more detailed discussions of these functions, and their relations to Goncharov's 
multiple polylogarithm, see Appendix B. It is clear from the definition that these functions 
form a shuffle algebra 

M{wi)M{w2) = Yl (^•^^) 

Using these functions we can easily integrate X^^^ and jj^^ . We illustrate this procedure 
explicitly for the integral 

'^'^ lnfi;f— -1^ +1^ ln(l--i;). (6.33) 



(6.34) 



First we can express all logarithms in terms of the G-functions we defined: 
In (v - 1^ + 1^ ln(l -v) = In ^1 - ^ ln(l - v) 

= G(A3,l;i;)+G(l,A3;i;). 

Then we get 

= / di; + / dv (6.35) 

.70 - Ai Jo u - Ai ^ ' 

= G(Ai,A3,l;l)+G(Ai,l,A3;l) 
= M(Ai,A3,l) + M(Ai,l,A3). 

All other integrals can be performed in exactly the same way, and we can hence express 
Xq^^ as a combination of M-functions. We find 

4'^(xi,X2)= (6.36) 

^|(iln2x2 + 4)M(Ai) + (-iln2x2-^)M(A2) -lnx2M(Ai,0) - 

lnx2M(Ai,l) +lnx2M(Ai,A3) + lna;2M(A2, O) +lnx2M(A2,l) - ln.X2M(A2, A3) + 
M(Ai,0,0) +M(Ai,0,l) -M(Ai,0,A3) +M(Ai,l,0) +M(Ai,l,l) -M(Ai,l,A3) - 
M(Ai,A3,0) -M(Ai,A3,l) +M(Ai,A3,A3) -M(A2,0,0) - M(A2, 0, l) + M(A2, 0, A3) - 

M(A2,1,0) -M(A2,1,1) +M(A2,1,A3) +M(A2,A3,0) + M(A2, A3, l) - M(A2, A3, A3) [. 
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Note that this expression is of uniform weight 3, as expected. 

The integration of X^^'* can be done in a similar way as for \ However, there is a 
slight complication. The function i^^^ contains polylogarithms of the form 

In order to perform the integration in terms of G-functions, we have to express these 
functions in terms of objects of the form G{. . . ; v). In Appendix G we show that the 
following identities hold: 

Li2 (^%;#^) = (6.38) 



v{v—l) 

-iln^xi + lna;2 lna;i - In^ X2 - G(0,0;?;) - G(0, 1; i;) + G (0, Ai;t;) + G(0, X2;v) - 

0; v) - \-v)+G (1, \i;v)+G (1, A2; v) + G (A3, i);v) + G (A3, 1; v) - 
G (A3, Xi;v)-G (A3, A2; + G(0; v) \nx2 + G(l; v) \nx2-G (A3; v) \nx2-M (0, Ai) - 
M (0, A2) + M (Ai, 1) + M (A2, 1) - M (A3, 0) - M (A3, 1) + M (A3, Ai) + M (A3, A2) - ^. 

Lis = (6.39) 

- In^ xi - ^G(0; v) In^ xi - iG(l; v) In^ xi + (A3; v) In^ xi - \M (A3) In^ xi + 

G(0; ?;) In X2 Inxi + G(l; v) lnx2 Inxi — G (A3; ?;) Inx2 Inxi + lna;2M (A3) Inxi + ^tt^ Inxi — 

G(0; ^;) In^ X2 - G(l; v) In^ ^2 + G (A3; v) In^ ^2 - Itt'^G{Q- v) - 1'it'^G{1- v) + ^tt^G (A3; ^;) - 

G(0, 0, 0; z;) - G(0, 0, 1; t;) + G (0, 0, Ai; i;) + G (0, 0, A2; v) - G(0, 1, 0; v) - G(0, 1, 1; v) + 

G (0, 1, Ai; t-) + G (0, 1, A2; v) + G (0, A3, 0; f) + G (0, A3, 1; f ) - G (0, A3, Ai; i;) - 

G(0,A3,A2;i;) -G(l,0,0;i') -G(l,0,l;z;) +G(l,0,Ai;t;) + G(1,0,A2;«) - G(l,l,0;t;) - 

G(l,l,l;i;) + G(l,l,Ai;i;) + G(l,l,A2;^;) + G(l,A3,0;t;)+G(l,A3,l;^;)-G(l,A3,Ai;^;)- 

G (1, A3, A2; v) + G (A3, 0, i)-v)+G (A3, 0, 1; t;) - G (A3, 0, Ai; i;) - G (A3, 0, M^v) + 

G (A3, 1, 0; ^;) + G (A3, 1, V^v) - G (A3, 1, Ai; i;) - G (A3, 1, A2; i;) - G (A3, A3, 0; v) - 

G (As, A3, 1; i;) + G (A3, A3, Ai; t;) + G (A3, A3, A2; w) + G(0, 0; v) In X2 + G(0, l-v) In X2 - 

G(0, A3;i')lnx2 + G(l, 0; t;) In X2 + G(l, l;t;)lna:2 - G(l, Ag;!;) \nx2 - G{\?„Q-v) \nx2 - 

G{\?„l:v) \nx2 + G{X?,,\^;v)\nx2 - X2M {X?,) - ivr^M (A3) - G(0; z;)M (0, Ai) - 

G(l; v)M (0, Ai) + G (A3; t;) M (0, Ai) - M (A3) M (0, Ai) - G(0; f )M (0, A2) - 

G(l; v)M (0, A2) + G (A3; v) M (0, A2) - M (A3) M (0, A2) + lna;2M (0, A3) + 

G(0; v)M (Ai, 1) + G(l; v)M (Ai, 1) - G (A3; v) M (Ai, 1) + M (A3) M (Ai, 1) + 

G(0; v)M (A2, 1) + G(l; v)M (A2, 1) - G (A3; t;) M (A2, 1) + M (A3) M (A2, 1) - 

G(0; i;)M (A3, 0) - G(l; v)M (A3, 0) + G (A3; t;) M (A3, 0) + lnx2M (A3, 0) - 

M (A3) M (A3, 0) - G(0; t;)M (A3, 1) - G(l; v)M (A3, 1) + G (A3; v) M (A3, 1) - 

M (A3) M (A3, 1) + G(0; t;)M (A3, Ai) + G{l-v)M (A3, Ai) - G (A3; ^;) M (A3, Ai) + 

M (A3) M (A3, Ai) + G(0; i;)M (A3, A2) + G(l; v)M (A3, A2) - G (A3; v) M (A3, A2) + 

M (A3) M (A3, A2) - lnx2M (A3, A3) - M (0, 0, Ai) - M (0, 0, A2) + M (0, Ai, 1) + 

M (0, A2, 1) - M (0, A3, 0) - M (0, A3, 1) + M (0, A3, Ai) + M (0, A3, A2) - M (Ai, 1, 1) - 

M (A2, 1, 1) - M (A3, 0, 0) + M (A3, 0, Ai) + M (A3, 0, A2) + M (A3, 1, 1) - M (A3, Ai, 1) - 

M (A3, A2, 1) + M (A3, A3, 0) + M (A3, A3, 1) - M (A3, A3, Ai) - M (A3, A3, A2) ) . 
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Using these identities we can express i^^^ completely in terms of G and M-functions, and 
perform the integration in exactly the same way as for Iq^ . The result is 

l['\x,,X2)= (6.40) 
^|(ln2x2 + ^)M(Ai,0) -|7r2M(Ai,l) + (i In^ X2 + ^)m(Ai, Ai) + 

i ln2 X2 + ^) M ( Ai , A2) + ( - ^ ln2 xi + In X2 In xi - 2 In^ a;2 - ^ ) M (Ai , A3) + 

- In^ X2 - ^ ) M(A2, 0) + |7r2M(A2, l) + ( - ^ In^ X2 - 4) ^(^2, Ai) + 

- i In^ X2 - ^)m(A2, A2) + (i In^ xi - \nx2 Inxi + 2 In^ X2 + ^) M(A2, A3) + 
-iln^xi +lnx2lnxi - In^ X2 - (A3, Ai) + 

iln^xi -lnx2lnxi + In^ X2 + ^) M(A3, A2) - 21nx2M(0, Ai, Ai) + lnx2M(0, Ai, A3) + 

21nx2M(0,A2,A2) - lnx2M(0, A2, A3) + lnx2M(0, A3, Ai) - lnx2M(0, A3, A2) - 
21nx2M(Ai,0,0) - lnx2M(Ai, 0, Ai) - lnx2M(Ai, 0, A2) + 21nx2M(Ai, 0, A3) + 
21nx2M(Ai,l,l) +lnx2M(Ai,l,Ai) - lnx2M(Ai, 1, A2) - lnx2M(Ai, 1, A3) - 
lnx2M(Ai,Ai,0) +lnx2M(Ai,Ai,l) + lnx2M(Ai, Ai, A3) - lnx2M(Ai, A2, 0) - 
lnx2M(Ai,A2,l) +lnx2M(Ai,A2,A3) + 21nx2M(Ai, A3, O) - lnx2M(Ai, A3, l) + 
lnx2Af(Ai,A3,Ai) + lnx2M(Ai, A3, A2) - 21nx2M(Ai, A3, A3) + 2 In X2M(A2, 0, O) + 
lnx2M(A2,0,Ai) +lnx2M(A2,0,A2) - 2 lnx2M(A2, 0, A3) - 2 In X2M(A2, 1, l) + 
lnx2M(A2,l,Ai) -lnx2M(A2,l,A2) + lnx2M(A2, 1, A3) + lnx2M(A2, Ai, O) + 
lnx2M(A2,Ai,l) -lnx2M(A2,Ai,A3) + lnx2M(A2, A2, 0) - lnx2M(A2, A2, l) - 
lnx2M(A2,A2,A3) - 21nx2M(A2, A3, O) + lnx2M(A2, A3, l) - lnx2Af(A2, A3, Ai) - 
lnx2M(A2,A3,A2) + 21nx2M(A2, A3, A3) - lnx2M(A3, 1, Ai) + lnx2Af (A3, 1, A2) - 
lnx2M(A3,Ai,l) +21nx2M(A3,Ai,Ai) - lnx2M(A3, Ai, A3) + lnx2M(A3, A2, l) - 
21nx2M(A3,A2,A2) + lnx2M(A3, A2, A3) - lnx2M(A3, A3, Ai) + lnx2M(A3, A3, A2) - 
2M(0,0,Ai,Ai) +2M(0,0,A2,A2) -M(0,Ai,0,Ai) -M(0,Ai,0,A2) +M(0,Ai,l,Ai) - 
M(0,Ai,l,A2) +2M(0,Ai,Ai,l) -2M(0,Ai,Ai,A3) -M(0,Ai,A3,0) -M(0,Ai,A3,l) - 
M(0, Ai, A3, Ai) + M(0, Ai, A3, A2) + M(0, A2, 0, Ai) + M(0, A2, 0, A2) + M(0, A2, 1, Ai) - 
M(0, A2, 1, A2) - 2M{0, A2, A2, 1) + 2M(0, A2, A2, A3) + M{0, A2, A3, O) + M{0, A2, A3, l) - 
M(0, A2, A3, Ai) + M(0, A2, A3, A2) - M(0, A3, 0, Ai) + M(0, A3, 0, A2) - M(0, A3, 1, Ai) + 
M(0, A3, 1, A2) - M(0, A3, Ai, 0) - M{0, A3, Ai, l) + M(0, A3, A2, O) + M(0, A3, A2, l) + 
2M(Ai, 0, 0, 0) - M(Ai, 0, 0, Ai) - M(Ai, 0, 0, A2) - M(Ai, 0, 0, A3) - 2M(Ai, 0, 1, l) + 
M(Ai,0,1,A3) +M(Ai,0,Ai,l) -M(Ai,0,Ai,A3) +M(Ai,0,A2,l) -M(Ai,0,A2,A3) - 
2M(Ai,0,A3,0) - 2M(Ai, 1,0,1) + M(Ai, 1, 0, A3) - 2M(Ai, 1, 1, 0) - 4M (Ai, 1, 1, l) - 
M(Ai,l,l,Ai) +M(Ai,l,l,A2) +3M(Ai,l,l,A3) -M(Ai,l,Ai,l) +Af(Ai,l,Ai,A3) + 
Af (Ai, 1, A2, 1) - A^(Ai, 1, A2, A3) + Af (Ai, 1, A3, 0) + 3Af (Ai, 1, A3, l) + Af (Ai, 1, A3, Ai) - 
M{Xi, 1, A3, A2) - 2A4'(Ai, 1, A3, A3) + A^(Ai, Ai, 0, 0) + A^(Ai, Ai, 0, l) - A^(Ai, Ai, 0, A3) + 
M{Xi, Ai, 1, 0) - Af(Ai, Ai, 1, 1) + Af(Ai, Ai, 1, A3) - Af(Ai, Ai, A3, 0) + M{Xi, Ai, A3, l) + 
A^(Ai,Ai,A3,A3) +A^(Ai,A2,0,0) +A^(Ai,A2,0, l) - A^(Ai, A2, 0, A3) + A4-(Ai, A2, 1, 0) + 
A^(Ai, A2, 1, 1) - A^(Ai, A2, 1, A3) - A^(Ai, A2, A3, 0) - Af (Ai, A2, A3, l) +Af (Ai, A2, A3, A3) - 
3Af(Ai,A3,0,0) +3Af(Ai,A3,l,l) +Af(Ai,A3,l,Ai) - Af(Ai,A3,l,A2) - 
2A4'(Ai, A3, 1, A3) + mIXi, A3, Ai, 1) + M{Xi, A3, Ai, A3) - M{Xi, A3, A2, l) + 
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M(Ai, A3, A2, A3) - 2M(Ai, A3, A3, 1) + M(Ai, A3, A3, Ai) + M(Ai, A3, A3, A2) - 
2M(A2, 0,0,0) +M(A2,0,0,Ai) +M(A2, 0,0, A2) +M(A2, 0,0, A3) +2M(A2, 0,1,1) - 
M(A2,0,1,A3) -M(A2,0,Ai,l) +M(A2,0,Ai,A3) -M(A2,0,A2,1) +M(A2,0,A2,A3) + 
2M(A2, 0, A3, 0) + 2M(A2, 1, 0, l) - M(A2, 1, 0, A3) + 2M(A2, 1, 1, 0) + 4M(A2, 1, 1, l) - 
M(A2, 1, 1, Ai) + M(A2, 1, 1, A2) - 3M(A2, 1, 1, A3) - M(A2, 1, Ai, l) + Af (A2, 1, Ai, A3) + 
M(A2, 1, A2, 1) - M(A2, 1, A2, A3) - M(A2, 1, A3, 0) - 3M(A2, 1, A3, l) + M(A2, 1, A3, Ai) - 
M(A2, 1, A3, A2) + 2M(A2, 1, A3, A3) - M(A2, Ai, 0, 0) - M(A2, Ai, 0, l) + M(A2, Ai, 0, A3) - 
M(A2,Ai,l,0) -M(A2,Ai,l,l) +M(A2,Ai,l,A3) +M(A2,Ai,A3,0) +M(A2,Ai,A3,l) - 
M(A2,Ai,A3,A3) -M(A2,A2,0,0) -M(A2,A2,0,1) +M(A2,A2,0,A3) -M(A2,A2,1,0) + 
M(A2, A2, 1, 1) - M(A2, A2, 1, A3) +M(A2, A2, A3, 0) -M(A2, A2, A3, l) -M(A2, A2, A3, A3) + 
3M(A2,A3,0,0) -3Af(A2,A3,l,l) +M(A2,A3,l,Ai) -M(A2,A3,1,A2) + 
2M(A2, A3, 1, A3) + M(A2, A3, Ai, 1) - M(A2, A3, Ai, A3) - M(A2, A3, A2, l) - 
M(A2, A3, A2, A3) + 2M(A2, A3, A3, 1) - M (A2, A3, A3, Ai) - M(A2, A3, A3, A2) - 
M(A3,0,0,Ai) +M(A3,0,0,A2) -M(A3,0,Ai,0) -M(A3,0,Ai,A3) +M(A3,0,A2,0) + 
M(A3, 0, A2, A3) - M(A3, 0, As, Ai) + M(A3, 0, A3, A2) + M(A3, 1, 1, Ai) - M(A3, 1, 1, A2) + 
M(A3, 1, Ai, 1) - M(A3, 1, Ai, A3) - M(A3, 1, A2, l) + M(A3, 1, A2, A3) - M(A3, 1, A3, Ai) + 
M(A3,1,A3,A2) -M(A3,Ai,0,0) -M(A3,Ai,0,A3) +M(A3,Ai,l,l) -M(A3,Ai,l,A3) + 
2M(A3,Ai,Ai,A3) -M(A3,Ai,A3,0) -M(A3,Ai,A3,l) +2M(A3,Ai,A3,Ai) + 
M(A3, A2, 0, 0) +M(A3, A2, 0, A3) - M(A3, A2, 1, 1) +M(A3, A2, 1, A3) - 2M(A3, A2, A2, A3) + 
M(A3, A2, A3, 0) + M(A3, A2, A3, 1) - 2M(A3, A2, A3, A2) - M(A3, A3, 0, Ai) + 
M(A3, A3, 0, A2) - M(A3, A3, 1, Ai) +M(A3, A3, 1, A2) - M(A3, A3, Ai, O) -M(A3, A3, Ai, l) + 
2M(A3,A3,Ai,Ai) +M(A3,A3,A2,0) +M(A3,A3,A2,1) -2M(A3,A3,A2,A2) + 
M[X2)(^- lln^ xi +5lna;2ln^xi -ln^X2lnxi -gTT^lnxi + In^ 0:2 + ^tt^ lnx2 - Ca) + 



M 



i)^gln^xi — ilnx2ln^a;i + 111^X2 In xi +i7r^lnxi — In^ X2 — ^tt^ lnx2 + Ca^ • 



Note that is of uniform weight 4, as expected. 



6.2 Evaluation of the Mellin-Barnes integral in Region II 

In this section we evaluate the pentagon in Region II. Since the Regions 11(a) and 11(b) 
are related simply by xi X2, we only concentrate on Region 11(a). The procedure is 
very similar to Region I, i.e. we start by deriving a twofold Mellin-Barnes representation 
for the pentagon in this region, which we then reduce to an Euler-type integral, and finally 
we express the result in terms of Goncharov's multiple polylogarithms. Performing this 
rescaling (4.5) in the Mellin-Barnes representation (6.2) and after the change of variable 



Z3 Z — Zi — Z2, 

(6.41) 

Z4 — >■ —z — 2z' + zi, 
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we find 

7f(l,l,l,l,l;Qf)= ^ 



r(l - 2e) 



(27ri)4 

X r i-zi) r (-Z2) r (zi + ^2 + 1) r (-e - z' - 1) r (e - z + zi - / + 2) 
X r (-Z - Z2 - ^' + 1) r (zi + Z2 - z') r (-e + z + z' - 1) r (-zi + z' + 1) 

X r (z + Zl + 2Z2 + 2 (-Z1 - Z2 + z')) . 

(6.42) 

We think of the integration over z as the last one, and we analyze how poles in r(. . .+z) with 
leading behavior A^^ might arise. There is only one possibility, coming from the product 
r (-e - z' - 1) r (-e + z + z' - 1). Taking the residues at z' = -1 - e + n', n' G N, we find 

A-^r (-e + z + z' - 1) ^ A-^r (-2e - 2 + z + n') . (6.43) 

Taking the residues at z = 2 + 2e — n — u', n G N, we find 

A^r (-2e - 2 + Z + n') ^ ^-2-2e+n+n'_ (g_44) 

Since we are only interested in the leading behavior in A~^, we only keep the terms in 
n = ra' = 0. Hence, we find a twofold Mellin-Barnes representation for the pentagon in 
multi-Regge kinematics, 

lf(l,l,l,l,l;Qf) = rre^-^^^xS"^(K,ii,t2) (6-45) 

St2 

with 

^M^^^'^''^''^) = r(l + emi-e)2 (2;^ y_ dzidz2yry2^^r(-e-zi)r(-zi)^ 

X r (zi + 1) r (-e - Z2) r (-Z2) r (zi + Z2 + 1) r (e + zi + Z2 + 1) , 

(6.46) 

where yi and ?/2 are defined in Eq. (5.10). We checked that if we close the integration con- 
tours to the right, and take residues, we reproduce exactly the expression of the pentagon 
obtained from NDIM, Eq. (5.13). 

We now evaluate the Mellin-Barnes representation (6.45) and we derive an Eulcr in- 
tegral representation for the pentagon in multi-Regge kinematics. Let us concentrate only 
on the Mellin-Barnes integral. We can now use the identity (6.9), and we checked again 
numerically that we can exchange the Euler and the Mellin-Barnes integration. We find 

x{l-v{l- y2)y-' r (-e - Zl) r i-zi f F (zi + 1)^ . 
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Finally, we close the zi-contour to the right and take residues at zi = ni — e, ni G N. As in 
the case of 
series in e, 



the case of X^g, we can sum up the series of residues and expand the integrand in a power 



We find 



and 



ti,t2) = li''"-\yi,y2) + e y2) + 0{e^). (6.48) 

Jo \yiv^ - yiv -y2V + v-l) 

4"'\vum) = /' iv , , ^-'"fa-^-") (6.50) 

where j^^^ and j^^^ are functions depending on (poly) logarithms of weight 2 and 3 respec- 
tively in yi, ?/2 and f (See Appendix F for the explicit expressions). Note that this implies 
that XQ^^"^(yi, 2/2) and l[^^°'\yi,y2) will have uniform weight 3 and 4 respectively, as ex- 
pected. Furthermore note that the poles in t; = and v = 1 have cancelled out. We 
need however still to be careful with the quadratic polynomial in the denominator of the 
integrand, since it might vanish in the integration region. We will analyze this situation in 
the rest of this section. 

We know already that the phase space boundaries in Region 11(a) require 

-^+^>1. (6.51) 

Since this implies X2 > 1 and X2 > x\, wc get from Eq. (5.10) that < 2/1,2/2 < 1- Wc now 
turn to the quadratic denominator in Eqs. (6.49) and (6.50). The roots of this quadratic 
polynomial are 



(6.52) 



A'l = A'i(2/i, y2) = ^ (^-1 + yi + y2 - \[>C^ 



A2 = Ai(yi, ^2) = ^ ^-1 + 2/1 + 2/2 + 



where denotes the Kallen function 

A'k = ^'K{yi,y2) = A(-2/i,2/2, 1) = l + yf + y2 + 2yi- 2y2 + 2yiy2. (6.53) 

First, let us note that A'^ (2/1,1/2) > everywhere in Region 11(a), and hence the square 
root in Eq. (6.52) is well-defined. Second, it is easy to show that we have, 

-1< a; (2/1,^2) <0 and l<A'2(2/i,y2)<2. (6.54) 

For later convenience, let us note at this point the following useful identities 

-1 



A'iA2 
A'i + A'2 

A'l - A2 



y2 

-1 + yi - y2 



yi 



(6.55) 



yi 
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From Eq. (6.54) it follows now immediately that the quadratic denominators in Eqs. (6.49) 
and (6.50) do not vanish in the whole integration v G [0, 1], and hence all the integrals in 
Eqs. (6.49) and (6.50) are convergent. Using partial fractioning and the relations (6.55) we 
can write 

'{yi,y2) =-7= dv — - — + -?^ di; — — — , (6.56) 



and 



^r^(yi,y2) =^f d.^^^^i^if^ + ^ f d.^^^^i^if^. (6.57) 

V^- -^0 V-\2 y/X^ Jo v-Xi 

Let us conclude this section by introducing the function 

^ A^(yi,y2) = = . ) = A3. (6.58) 

1 — ^/2 i — tl/t2 

This function appears in j^^^ and j^^^ through the logarithm 

In (1 - t; (1 - y2)) = In (l - ^) = £ (6.59) 

We now evaluate the integrals I^^^"'^ and explicitly. For Iq^""^ , we find 

Xf'^^(yi,y2)= (6-60) 
-iln^yi -zj)M(A;) + (iln2yi + nl) M {\'^) - Iny^ M {X[,0) - 



Inyi M{X[, 1) + Inyi M{X[, A'3) + Inyi M(A'2,0) + Inyi M{X'^, l) - Inyi M{X'^, A'3) - 

m(a;,o,o) -m(a;,o,i) +M(A;,o,A^j) -m(a;,i,o) -m(a;,i,i) +m(a;,i,a;3) + 
m(a;,a^j,o) +m(a;,a'3,i) - m(a;,a;3,a^j) +m(a'2,o,o) +m(a'2,o,i) - m(a'2,o,a^j) + 

M(A'2, 1, 0) + M(A^, 1, 1) - M(A'2, 1, A0 - M(A^, A'3, 0) - M(A^, A'3, l) + M(A^, A^j, A'3) | . 

Note that this expression is of uniform weight 3, as expected. 

The integration of x|^^"^ can be done in a similar way as for Tq^"'\ However, there is 
again a slight complication. The function j^^^ contains poly logarithms of the form 

(y2v-v + l\ 
Lin —( TT — . (6.61) 

In order to perform the integration in terms of G-functions, wc have to express these 
functions in terms of objects of the form G{...\v). In Appendix H we show that the 
following identities hold: 

(ifHgr) = (6-62) 
-1 In^ (1^) - G(0, 0; v) - G(0, l;v)-G (0, A'^, I) + G (0, X'^-v) - G (0, A'2, 1) + 
G (0, A'2; v) - G(l, 0; v) - G(l, l]v)+G (1, A'^; v) + G (1, A'2; v) + G {X\, 1, 1) + G (A'2, 1, 1) - 
G (A'3, 0, l)+G (A'3, 0; v)-G (A'3, 1, 1) + G (A'3, l-v)+G (A'3, A'^, 1) - G (Aj,, X!^-v) + 
G (A^5, A'2, 1) - G (A^, A'2; v) - G(0; v) Inyi - G(l; v) Inyi + G (A'3; v) Inyi - Inyi lny2 - ^ • 
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(ff^) = (6.63) 
i In^ (f ) + iG(0, 1) In^ ) - iG(0; In^ (f ) - t;) In^ ) - 

iG (A'3, 1) In^ {^)+lG (A'3; In^ {^) + ^tt^ In + Itt^GCO, 1) - In^O; v) - 

^7r2G(l; t;) - ^tt^G (A'3, 1) + ^tt^G (A'3; ^) + G(0, 1)G (0, A'^, 1) - G(0; (0, A'^ , 1) - 

G{l;v)G (0, A'l, 1) - G (A'3, 1) G (0, A;, 1) + G {\i,;v) G (0, A'^, 1) + G(0, 1)G (0, A'2, 1) - 

G(0; ?;)G (0, A'2, 1) - G(l; 7;)G (0, A'2, 1) - G (A'3, 1) G (0, A'2, 1) + G (A'3; t;) G (0, A'2, 1) - 

G(0, 1)G (A;, 1, 1) + G(0; v)G {\[, 1, 1) + G(l; t;)G (A;, 1, 1) + G (A'3, 1) G (A'„ 1, 1) - 

G (A^; z;) G (A'^, 1, 1) - G(0, 1)G (A'2, 1, 1) + G(0; v)G (A'2, 1, 1) + G(l; ^;)G (A'2, 1, 1) + 

G (A^, 1) G (A'2, 1, 1) - G (A^; t;) G (A'2, 1, 1) + G(0, 1)G (A'3, 0, 1) - G(0; ^;)G (A'3, 0, 1) - 

G(l; z;)G (A'3, 0, 1) - G (A'3, 1) G (A'3, 0, 1) + G (A'3; t;) G (A'3, 0, 1) + G(0, 1)G (A;^, 1, 1) - 

G(0; 7;)G (A'3, 1, 1) - G(l; ^;)G (A'3, 1, 1) - G (AJ,, 1) G (A'3, 1, 1) + G (A;,; t;) G (A'3, 1, 1) - 

G(0, 1)G (A'3, a;, 1) + G(0; v)G (A'3, A'^, 1) + G(l; t;)G (A;,, A;, 1) + G (A'3, 1) G (A'3, A'^, 1) - 

G (A'3; G (A'3, a;, 1) - G(0, 1)G (A'3, A'2, 1) + G(0; .;)G (A'3, A'2, 1) + G(l; t;)G (A'3, A'2, 1) + 

G (A^5, 1) G (A'3, A'2, 1) - G (A'3; G (A^, A'2, 1) + G(0, 0, 0, 1) - G(0, 0, 0; - G(0, 0, 1; z;) - 

G (0, 0, A'l , 1) + G (0, 0, X[-v) - G (0, 0, A'2, 1) + G (0, 0, A'2; i;) - G(0, 1, 0; v) - G(0, 1, 1; v) + 

G (0, 1, X[;v) + G (0, 1, A'2; v) + G (0, A;, 1, 1) + G (0, A'2, 1, 1) - G (0, A'3, 0, 1) + 

G (0, A'3, 0; v)-G (0, A'3, 1, 1) + G (0, A'3, 1; z;) + G (0, A'3, A;, 1) - G (0, A'3, A;; t;) + 

G (0, A'3, A'2, 1) - G (0, A'3, A'2; t-) - G(l, 0, 0; t;) - G(l, 0, 1; t-) + G (1, 0, X[;v) + 

G (1, 0, A'2; t-) - G(l, 1,0; t-) - G(l, 1, 1; i;) + G (1, 1, X[;v) + G (1, 1, A'2; t-) + G (1, A^^, 0; v) + 

G (1, A'3, 1;^)-G(1, A'3, A',; i;)-G(l, A'3, A'2; ^;)-G(A'i, 1,1,1) -G (A'2, 1,1,1) - 

G (A^5, 0, 0, 1) + G (A'3, 0, 0-v)+G (A'3, 0, 1; t;) + G (A^^, 0, A'^, 1) - G (A'3, 0, A'^; ^;) + 

G (A'3, 0, A'2, 1) - G (A'3, 0, A'2; v) + G (A'3, 1, 0;v) + G (A'3, 1, 1, 1) + G (A:3, 1, 1; v) - 

G (A'3, 1, X[;v) - G (A'3, 1, A'2; v)-G (A'3, A;, 1, 1) - G (A'3, A'2, 1, 1) + G (A'3, A'3, 0, 1) - 

G (A'3, A'3, 0; v) + G (A'3, A'3, 1, 1) - G (A;,, A:3, 1; t;) - G (A:3, A'3, A'l, 1) + G (A:3, A:3, X[;v) - 

G (A^, A'3, A'2, 1) + G (A'3, A'3, A'2; t;) + G(0, 0, 1) In (y,) - G(0, 0; v) In (y,) - G(0, 1; v) In (yi) - 

G (0, A'3, 1) In (j/i) + G (0, A'3; v) In (yi) - G(l, 0; v) In (yi) - G(l, 1; v) In (yi) + 

G (1, A'3; ^;) In (t/i) - G (A^^, 0, 1) In (yi) + G (A'3, 0; ^;) In (yi) + G (A^, 1; v) In (yi) + 

G (A:3, a;,, 1) In (yi) - G (A'3, A'3; t;) In (yi) + G(0, 1) In (yi) In (y2) - G(0; v) In (yi) In (ys) - 

G{l;v) In (yi) In (y2) - G (A^, 1) In (yi) In (y2) + G (A'3; v) In (yi) In (y2) ■ 



Using these identities we can express j^^^ completely in terms of G and M- functions, and 
perform the integration in exactly the same way as for X^^^"^ . The result is 

li'"'\yi,y2)= (6-64) 

I ( - 2 ln2 yi - ln2 y2 - ^ ) M (A; , o) + ( - In^ yi - In^ ya + 4) M (A'^ , l) + 

(-iln^yi -^)M(A'„Ai) + (-iln2yi -^)m(A'i,A0 + 
(fln^yi +iln2y2 + ^)M(A'i,A'3) + (21n2yi + In^ y2 + ^) M(A'2, O) + 
(in^yi +ln2y2-^)M(A'2,l) + (^ln2yi + if ) M(A^, A'^) + In^ yi + ^ 

( - I ln2 yi - 1 In^ y2 - 4) M(A'2, A0 + ( - i In^ y, - i In^ y2 - ^) M(A'3 




-)m(a^,a0 

y2-^)M(A'3,A'0 + 
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(i ln2 yi + i ln2 + ^)m(A'3, - 21nyi Af(0, A^, A;) - Inyi M(0, A;, A'3) + 

21nyiM(0,A'2,A'2) +lnyiM(0,A'2,A;,) -lnyiM(0,A'3,A;) +lnyiM(0,A^5,A'2) - 
41nyiM(A;,0,0) -21nyiM(A;,0,l) -lnyiM(A;,0,A;) -lnyiM(A'i,0,A'2) + 
21nyiM(A'i,0,A'3) -21nyiM(Ai,l,0) +lnyiM(A'i,l,A'i) -lnyiM(A'i,l,A0 + 
lnyiM(A'i,l,A'3) -ln2/iM(A'i,Ai,0) + Inyi M(A'i, A'^, l) + Inyi M(A'i, A'^, AQ - 
ln2/iM(A'i,A^,0) -lnyiM(A'i,A^,l) + Inyi M(A'i, A'2, A^) + 2 Inyi M(Ai, A'3, 0) + 
In 2/1 M(A'i, A?5, 1) + Inyi M(A;, A'3, A^) + Inyi M(A;, A'3, A'2) - 21nyi M(A'i, A^,, A'3) + 
41nyiM(A'2,0,0) + 21nyi M(A'2,0, l) + Inyi M(A'2, 0, A;) + Inyi M(A'2, 0, A'2) - 
21nyiM(A'2,0,A'3) +21nyiM(A'2,l,0) +lnyiM(A'2,l,A'i) -ln|/iM(A'2,l,A'2) - 
Inyi M(A'2, 1, A'3) + Inyi M(A'2, A^, O) + Inyi M(A'2, A'^, l) - Inyi M(A'2, A'^, A0 + 
lnyiM(A'2,A'2,0) -ln2/iM(A'2,A^,l) -lnyiM(A'2,A'2,A^,) -21nyiM(A^,A'3,0) - 
Inyi M(A'2, A'3, l) - Inyi M(A'2, A;,, A'^) - Inyi M(A'2, A;,, A'2) + 21nyi M(A'2, Aj,, AQ - 
21nyiM(A^j,0,A'i) +21nyiM(A'3,0,A'2) -lnyiM(A^5,l,A'i) +lnyiM(A'3,l,A'2) - 
21nyiM(A;j,A;,0) - Inyi M(A'3, A^, l) + 21nyi M^A^j, A;, A;) + Inyi M(A;3, A;, A'3) + 
21nyi M(A^5, A'2, O) + Inyi M{\'^, A'2, l) - 21nyi M(A'3, A'2, A'2) - Inyi M(A'3, A'2, A^j) + 
lnyiM(A'3,A^j,A'i) -lnyiM(A(j,A'3,A'2) -2M(0,0,A;,A'i) +2M(0,0,A'2,A'2) - 
3M(0,Ai,0,Ai) -3M(0,Ai,0,A'2) -M(0,A'i,l,Ai) -3M(O,Ai,l,A0 -4M(0,A'i,A'i,0) - 
2M (0, A'l , a; , 1) + 2M (0, A'l , A'^ , A^j) - M (O, A; , A^, , O) - M (O, A'^ , A'3, l) + M (O, A; , A'3, A'/) + 
3M(0, A'l, A'3, A'2) +3M(0, A'2, 0, A'l) +3M(0, A'2, 0, A'2) +3M(0, A'2, 1, A'^) +M(0, A'2, 1, A'2) + 
4M(0,A'2,A'2,0) +2M(0,A'2,A'2,l) -2M(0,A'2,A'2,A'3) +M(0,A'2,A;3,0) +M(0,A'2,A'3,l) - 
3M(0,A'2,A'3,A'i) -M(0, A'2, A'3, A'2) - M(0, A'3, 0, A'^) + M(0, A'3, 0, A'2) - M(0, A^j, 1, A^) + 
M(0, A^j, 1, A'2) - M(0, A^j, A'l, 0) - M(0, A'3, A'^, l) + M(0, A'3, A^, O) + M(0, A'3, A^, l) - 
4M(A'i, 0,0,0) -2M(A'i, 0,0,1) - 3M(A'i,0,0, A'^) - 3M(Ai, 0,0, A^) + M(A'i, 0,0, A^) - 

2m(a;, 0,1,0) -m(a;,o,i,a;;5) -2m(a;,o,a;,o) + m(a;,o,a;,i) +m(a;,o,a;,a;3) - 

2M (a; , 0, A'2 , 0) + M ( a; , 0, A'2 , 1) + M (A'l , 0, A'2 , a:3) - 2M (a; , 0, a;, , 1) + 2M (a; , 0, A^j , A'l) + 

2m(a;, o,A^5, A'2) -2m(a;, 1,0,0) -m(a;,i,o,aQ +2m(a'i,i,i,i) +m(a'i,i,i,a;) + 

3M(A'i, 1, 1, A'2) - 3M(A'i, 1, 1, A'3) +2M(A'i, 1, A'^, O) +3M(a;, 1, A;, l) -M(A'i, 1, A'^, A'3) - 
2M(A'i, 1, A^, 0) + M(a;, 1, A'2, 1) + M(a;, 1, A'2, A'3) - M(A'i, 1, A^j, 0) - 3M(a;, 1, A'3, 1) + 
M(A'i, 1, A'3, Ai) -M(Ai, 1, A^, A'2) +2M(A'i, 1, A^, A^) -M(Ai, Ai,0, 0) +3M(A'i, Ai,0, l) + 

M ( a; , a; , 0, + 3 Af (a; , a; , i , o) + 5m (a; , a; , i , i) - m (a; , a; , i , a^j) + m ( a; , a; , a^j , o) - 
M(A'i, a;, A'3, 1) - m(a;, a;, a:3, A'3) - m(a;, A'2, 0, 0) - m(a;, A'2, 0, 1) + m(a;, A'2, 0, a^,) - 
M(A'i,A'2,i,o) -m(a;,a'2,i,i) +m(a;,a'2,i,a'3) + m(a;,a'2,a'3,o) +m(a;,a'2,a'3,i) - 

M(A'i, A'2, A'3, A'3) -M(A'i, A'3, 0, 0) -2M(A;, A^, 0, l) +2M(A'i, A'3, 0, A'i)+2M(a;, A'3, 0, A'2)- 
2M(A'i, A'3, 1, 0) -3M(A'i, A'3, 1, 1) +M(A'i, A^j, 1, A'^) -M(A'i, A'3, 1, A'2) +2M(A'i, A^,, 1, A^,) + 
2M(A'i,A'3,Ai,0) +M(A'i,A^,A'i,l) -M(Ai,A^,A'i,A^) +2M(A'i,A'3,A^,0) - 

M(A'i,A^,A'2,i) -m(a;,a'3,a'2,a^5) + 2M(A;,A'3,Af5,i) -m(a;,a?j,a'3,a;) - 

M(a;,A'3,A;3, A'2) +4M(A'2, 0,0,0) +2M(A'2,0,0, l) + 3M(A'2, 0, 0, A'^) + 3M(A'2, 0, 0, A'2) - 
M(A'2,0,0,A^j) +2M(A'2, 0,1,0) +M(A'2, 0,1, A'3) +2M(A'2,0,A'i,0) -M(A'2,0,A'i,l) - 
M(A'2, 0, A'l, A^j) + 2M(A^, 0, A^, O) - M(A'2, 0, A^, l) - M(A'2, 0, A'2, A^,) + 2M(A^, 0, A^,, l) - 
2M(A'2, 0, \'^, a;) - 2M(A^, 0, \'^, AQ + 2M(A^, 1, 0, 0) + M(A'2, 1, 0, A'3) - 2M(A'2, 1, 1, l) - 
3M(A'2,1,1,A;) -M(A'2,1,1,A'2) +3M(A'2,1,1,A^) +2M(A'2,1,A;,0) -M(A'2,1,A;,1) - 
M(A'2, 1, A'l, AQ - 2M(A^, 1, A^, 0) - 3M(A^, 1, A^, l) + M(A'2, 1, A^, AQ + M(A^, 1, A^, O) + 
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3M(A'2, 1, A'3, 1) +M(A'2, 1, A'3, a;) -M(A'2, 1, A'3, AQ - 2M(A'2, 1, A^,, A:3) +M(A'2, A;, 0, 0) + 

m(a'2,a;,o, 1) -m(a'2,a;,o,a:3) +m(a'2,a;,i,o) +m(a'2,a;,i,i) -m(a'2,a;,i,a'3) - 

M(A'2, a;, A'3, 0) -M(A'2, a;, A'3, 1) +M(A'2, X[, A'3) +M(A'2, A'2, 0, 0) -3M(A'2, A'2, 0, l) - 
M(A'2, A^, 0, A^) -3M(A^, A^, 1, 0) -5M(A'2, A^, 1, l) +M(A'2, A^, 1, A^) -M(A^, A'2, A^, O) + 
M(A'2, A'2, A'3, 1) +M(A'2, A'2, A'3, A0 +M(A'2, A^,, 0, 0) +2M(A'2, A'3, 0, l) -2M(A'2, A^,, 0, A'^) - 
2M(A'2, A'3, 0, A'2) +2M(A'2, A^j, 1, 0) +3M(A'2, A'3, 1, l) +M(A'2, A'3, 1, A'^) -M(A'2, A^, 1, A0 - 
2M(A'2, A:3, 1, AQ - 2M(A'2, A:3, A;, O) + M(A'2, A'3, A;, l) + M(A'2, A'3, A;, A^) - 
2M(A'2, A'3, A'2, 0) - M(A'2, A^j, A'2, 1) + M(A'2, A;,, A'2, A;,) - 2M(A'2, A;3, A'3, l) + 
M(A'2, X's, A'3, A'/) +M(A'2, A'3, X's, A'2) - M(A'3, o', 0, A'^) +M(A^j, 0, 0, A'2) -3M(A'3, 0, A'^, O) - 
2M(A'3, 0, a;, 1) +M(A'3, 0, A'l, A^) +3M(A^, 0, A'2, 0) +2M(A^, 0, A'2, l) - M(A^„ 0, A^, A'3) - 

m(a;3,o,a;3,a;) + m(a;;5,o,a;^,a'2) + Af(A'3,i,i,A;) -m(a^,i,i,a'2) -2Af(A;^,i,A;,o) - 
m(a'3,i,a;,i) +m(a:3,'la;^a'3) +2M(A'3,i,A'2,o) +m(a:3,i,a'2,i) -m(a;5,'i,a'2,a:3) - 
m(a;3, 1, AJs, a;) +M(A'3, 1, A'3, A'2) -5m(a;3, a;, o, o) -4m(a'3, a;, 0, i) +2m(a;3, a;, o, a;) + 

2M(A'3, A'l, 0, A'2)+M(A^5, A'^, 0, A^) -4M{X'„ X[, 1, 0) -3M(A:3, A^, 1, l) +2M(A^j, A'^, 1, X[) + 

2M(A'3,A'i,1,A'2) +m(a'3,a;,i,a0 +4M(A'3,A;,A'i,0) +4M(A'3,A;,A'i,1) - 

2M(A'3,A'i,Ai,A'3) +M(A^,A'i,A^,0) +M(A'3,A'i,A^,l) -2M(A'3,A'i,A^,A'2) + 

5M(A'3, A'2, 0, 0) + 4M(A^5, A'2, 0, l) - 2M(A'3, A'2, 0, A;) - 2M(A'3, A'2, 0, A'2) - 

M(A'3, A'2, 0, A;3)+4M(a;;5, A'2, 1, 0)+3M(a;3, A'2, 1, 1) -2M(A^j, A'2, 1, A;)-2M(A^j, A'2, 1, A'2)- 

M(A:3, A2, 1, Xi,) - 4M(A'3, A'2, A'2, 0) - 4M(a:3, A'2, A'2, l) + 2M(A'3, A'2, A'2, A^,) - 

M(A'3, A'2, A:3, 0) - M(A'3, A'2, A'3, 1) + 2M(A'3, A'2, a:3, A^) - M(A'3, A^^, 0, A'^) + 

M(A'3, x'^, 0, A'2) -M(A^j, A^5, 1, a;) +m(a^5, x'^, 1, A'2) - m(a^j, A'3, a;, o) -m(a'3, A'3, a;, 1) + 

2M(A'3,A'3,Ai,A'i) +M(A^,A'3,A^,0) +M(A'3,A'3,A^,1) -2M(A'3,A'3,A^,A'2) + 

M(A'i) ( - I In^ yi + i Inys In^ Vi - In^ y2 Inyi - ^tt^ Inyi + i In^ y2 + ^tt^ In 7/2 - Cs) + 

M (A'2) (1 In^ yi - i In y2 In^ + In^ j/2 In j/i + ivr^ In yi - i In^ y2 - In y^ + Ca) | ■ 

Note that X^^'^^ is of uniform weight 4, as expected. 

7. Analytic continuation to the physical region 

The result in the physical region where all s-typc invariants are positive is obtained by 
performing analytic continuation on all s-type invariants according to the prescription 

(-s)^e-^'^s, (-si)^e-^'^si, (-S2) ^ e"^'^ S2. (7.1) 

The prescription for the transverse scale k is then fixed by Eq. (4.7) to be [38] 

(-k) ^ e"'^ At. (7.2) 

7.1 Analytic continuation of Region 11(a) 

In Region 11(a) the pentagon is expressed in terms of the dimensionless quantities y\ and 7/2 • 
The analytic continuation of y\ and 2/2 follows then directly from Eq. (5.10) and Eq. (7.2), 

(-yi) e"*''^! and y2 ^ y2- (7.3) 
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In the physical region, the pentagon can then be written as 



4S?(^>^i>^2,ii,i2) =rre^-^^4i;?(/.,ii,i2), (7.4) 



with 



T^l^^{^MM) = e^^^l\^^^\e-^-{-^)M,t2), (7.5) 



and 1^^^^ is given in Eq. (5.13) to all orders in e in terms of Appell and Kampe de 
Feriet functions. As we face the problem how of to perform the analytic continuation of 
these functions under the prescription (7.3), we argue in the following that the analytic 
continuation of the generalized hyper geometric functions that appear in l^^^i^ is trivial, 
i.e. , the hypergeometric functions stay real in the physical region. 

Indeed all the hypergeometric functions entering l^ii^ can be written as a nested 
sum, the inner sum being a one-dimensional hypergeometric function of 2F1 or 3F2 type 
depending on (— yi). For example, the first term in Eq. (5.13) can be written as 



^ y^' r(l - 2e) r(l + ef F4 (l - 2e, 1 - e, 1 - e, 1 - e; -yi, 
= - 1 yr r(l - 2e) r(l + ef ± ' e).. yj ^. \ 

^ n=0 ^ ^'"^ 

(7.6) 



Since we do not perform any analytic continuation in y2, the only way an imaginary part 
could arise is when we cross a branch cut during the analytic continuation in yi . Hence the 
analytic properties of Eq. (7.6) are determined by the cuts in the function 2F1 ^1 — 2e, 1 — 

e, 1 — e; — yi^ . We know that the 2F1 function has a branch cut ranging from 1 to +00. The 
convergence criterion for the i<4 function requires |yi| < 1, so we do not cross the branch 
cut when continuing yi along a half circle through the upper half plane, and so we do not 
change the Riemann sheet during the analytic continuation. Since the Euclidean region 
corresponds to the Riemann sheet where the hypergeometric function is real for (— yi) < 1, 
we conclude that the 2F1 function in Eq. (7.6) is real in the physical region, and so is the 
Appell function on the left-hand side of Eq. (7.6). A similar reasoning can be made for all 
other hypergeometric functions in Eq. (5.13), and so we can immediately write down the 
analytic continuation of Eq. (5.13) to the physical region where all s-type invariants are 
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positive, 



- ^ 2/^ e'"' r(l - 2e) r(l + ef F4 (l - 2e, 1 - e, 1 - e, 1 - e; -^1,^2) 
+ jL e^T6 r(i + e) r(l - e) F4 1 - e, 1 - e, 1 + e; -yi, ys) 



^2,1 



l + (5 l + (5-e 



d5 



l + 1-e l + e + 5 - 



|<5=0 



+ [M-yi) - iT^ + ip{l - e) - V'(-e)] -^4(1,1 - e, 1 + e, 1 - e--yi,y2) 



9 p2,l 



l + 1 + e l + e + (5 - 



- yi,y2 



1 5=0 



+ [In(-yi) -m + ^'(1 + e) - V;(-e)] ^4(1, 1 + e, 1 + e, 1 + e; -yi, 2/2) 



(7.7) 



We checked explicitly that all the poles in e cancel out when expanding ^Tp^y""* as a Laurent 
series in e. Note that in Eq. (7.7) the same hyper geometric functions appear as in the 
Euclidean region, Eq. (5.13). We can therefore easily expand Eq. (7.7) in e in exactly 
the way as we did in the Euclidean region, and we will obtain a Laurent scries whose 
coefficients are combinations of real A^-functions. The imaginary parts arise order by 
order in the Laurent series only through the explicit dependence on m in Eq. (7.7). Note 
that if we had worked with the expression of the pentagon in terms of Goncharov's multiple 
poly logarithms, Eqs. (6.60) and (6.64), the complicated analytic structure of these functions 
would imply many spurious imaginary parts arising from individual polylogarithms, and 
all the spurious imaginary parts would need to cancel out in the final answer. It is therefore 
more natural to express the result in the physical region in terms of M. functions rather 
than the Goncharov polylogarithms. 

7.2 Analytic continuation of Region I 

In Region I the pentagon is expressed as a function of xi and X2, whose analytic continu- 
ation follows directly from Eq. (4.6) and Eq. (7.2), 

(— xi) — e+^'^xi and X2 — 6+'^^ X2, 

and the pentagon in Region I can be expressed as 

-^Ss(*'-5i'*2,ii,t2) = rre^^' 

with 



S1S2 



e-^xiS(e— (-«:),ti,t2), 



(7.8) 

(7.9) 
(7.10) 



-37- 



and 1^1) is given in Eq. (5.29). Using the same argument as in Region 11(a), we see that the 
hypergeometric functions stay real even in the physical region, because the convergence of 
the hypergeometric functions requires \xi\ < 1, i = 1,2. Alternatively, we could also start 
from the expression of the pentagon in the physical Region 11(a), Eq. (7.7), and perform 
the analytic continuation to Region I following the techniques described in Appendix E. 
We find that the two results are consistent. The final expression of the pentagon in the 
physical region I is then 

- ^ e-'^^ {-xi)~^ {-X2)~^ r(l - 2e) T{1 + ef F^i} - 2e, 1 - e, 1 - e, 1 - e; -xi, -X2) 
+ ^ r(l - e)r(l + e) e'^' ^4(1, 1 + e, 1 + e, 1 + e; -xi, -X2) 



1-e l + e + 5 -1 + 5 



Xi,-X2 



\s=o 



+ [ln(-X2) + ivr + - e) - V'(-e)] ^4(1, 1 - e, 1 - e, 1 + e; -xi, -X2) 



_J_ ti 0,3 [ - - 
^ l-et2 I 2 2 - 6 



1 1 1 1 l-€ 1-e 



Xl 
X2 



1 ( ,^-J d 2,1 l + 5l-e + 5 



- - 1 
1+e l-e+5 -1+5 



- Xi,-X2 



I (5=0 



+ [ln(-X2) + m + tpil - e) - V'(e)] -^4(1, 1 - e, 1 + e, 1 - e; -xi, -X2) 



e t 



1^0,3 



1 + e i2 



2,0 



2 2 + e 



1 1 1 1 1-e 1 + e 



- xi, 




(7.11) 



We checked explicitly that all the poles in e cancel out when expanding I^^^^ as a Laurent 
series in e. 



8. Conclusions 



In this paper we computed the one-loop scalar massless pentagon integral in D = 6 — 2e 
dimensions, /f "^*^, in the limit of multi-Regge kinematics. This integral first contributes to 
the parity-odd part of the one-loop J\f = 4 five-point MHV amplitude at 0{e) - see (1.3). 

It is well known that loop integrals are intimately related to multiple generalised hy- 
pergeometric series and generalised polylogarithms. We have exploited two completely 
different approaches - Negative Dimensions and Mellin-Barnes integration - to extract the 
leading two terms of the Laurent expansion as e — for the in the high energy limit 

defined by, 

S > Si, S2 > -ti,-t2. (8.1) 
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In this limit, the pentagon integral reduces to double sums or equivalently two- fold Mellin- 
Barnes integrals. Equations (5.21), (5.22), (5.23) and (6.45), (6.48), (6.60), (6.64) give 
explicit expressions for the 0{€^) and 0{e) contributions to in terms of M. functions 

and Goncharov's multiple polylogarithms for kinematic region 11(a) when (— ti) < (— i2) 
and — a/— sti/siS2 + \/ —st2/ siS2 > 1. Analogous expressions in region I defined by 
■\/ —sti/ siS2 + \/—st2/siS2 < 1 can be derived through the analytic continuation of Sec- 
tion E.3 and are given in Eqs. (6.6), (6.12), respectively. 

By determining the 0(e) contribution to one therefore gains knowledge of m^^^ 

through to 0{e^) which is necessary for studies of the iterative structure of A/' = 4 SYM 
amplitudes beyond one-loop. One immediate application is the extraction of the one-loop 
gluon-prodTiction vertex through to C(e^) and the iterative construction of the two-loop 
gluon-production vertex through to finite terms which is described in a companion paper [1]. 

Finally, we note that one- loop polygons are present in the 0{e) contributions [22] 
to one-loop = 4 MHV amplitudes. In particular, the D = 6 — 2e hexagon integral 
contributes to the parity-even part [11, 22] of the one-loop six-point MHV amplitude, 
while D = 6 — 2e pentagon integrals contribute to the parity-odd part [12, 22]. It seems 

(2) 

reasonable to conjecture that the analytic form of the remainder function may be 
somewhat linked to the properties of the special functions occurring in the hexagon one- 
loop integral. The analytic methods we have used for evaluating the pentagon integral 
in the high energy limit may also be applied to the hexagon integral and may ultimately 

(2) 

provide information on the form of the i?g remainder function. 
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A. Nested hcirnionic sums 

The nested harmonic sums have been defined in Eq. (5.18). The S and Z sums fulfill an 
algebra. Let us illustrate this with a simple example: 

n n ^ 

S^{n)Sj{n) = E E 7-7J 

jtl = l k2 = l '^l'^2 
n ^ ki ^ n ^ k2 -. rt _ 

Ei ^ — ^ i ^ — ^ 1 ^ — ^ i ^ — ^ 1 

fci = l 1 k2 = l ^2 k2 = l ^2 ki=l 1 k=l 

= Sij{n) + Sji(n) - Si+j{n). 



(A.l) 
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A similar result can be obtained for the Z-sums^, 

Zi{n)Zj{n) = Zij{n) + Zji{n) + Zi+j{n). (A.2) 

For sums of higher weight, a recursive application of the above procedure leads then to 
the reduction of any product of 5 or Z sums to a linear combination of those sums. 
Furthermore, S and Z sums can be interchanged, e.g. 



k=i k=i e=i k=i \ e=i i 

k=l k=l 



(A.3) 



For n — GO, the Euler-Zagier sums converge to multiple zeta values, 

lim Zmi...mkin) =Cimk,...,mi). (A.4) 

n — >oo 

In Ref. [36] Moch et al. introduced generalizations of the S and Z-sums to make them 
dependent on some variables, 



Si{n]x) = Zi{n;x) = ^ 



k 



k=l 

x" 



Sij{n;xi,...,xe) = ^ '-^ Sj{k;x2, ■ ■ ■ ,Xn), (A.5) 



n I, 



Zij{n; xi,...,Xi) = ^ ^ Zj{k - 1; X2, • • • , Xn). 
k=l 

Those sums share of course all the properties of the corresponding number sums introduced 
in the previous paragraph. In particular, for n — ^ oo, the Z sums converge to Goncharov's 
multiple poly logarithm, 

lim Zmi,„rnk{n;Xi, . . . ,Xk) ='Limk,...,mi{xi, . . . ,Xk). (A.6) 

n — >oo 

These multiple polylogarithms will be reviewed in the next section. 

B. Goncharov's multiple polylogarithm 

Let us define [36] 

n{a) = (B.l) 

and iterated integrations by 

/•^ di /"* 

/ n{ai) o . . . o n{an) = / n{a2)o ...on{an). (B.2) 

Jo Jo t~ fli Jo 

®Note the sign difference with respect to Eq. (A.l). 
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We can then define the G and M-functions as 

G{ai,. . . ,an;z) = I 0(ai) o . . . oO(a„), 

(B.3) 

M(ai, . . . ,a„) =G(ai, . . . ,a„; 1) = / 0(ai) o . . . o 0(a„). 

Jo 

The M-functions are nothing but Goncharov's multiple polylogarithm, up to a sign^, 

M(ai,...,an) = {-lfUmk^,„^mi{zk,---,Zl), (B.4) 

where k is the number of nonzero elements in (oi, . . . , a„). 

Iterated integrals form a shuffle algebra, and hence wc can immediately write 



G{wi\ z) G{w'2\ z) = ^ G{w\z), 
M{wi)M{w2) = ^ M{w). 



(B.5) 



It is easy to see that in general* 



or equivalently 



lim M(ai, . . . , a„) = GO, (B.6) 

ai— >1 



lim G(l,a2, . . . ,a„;z) = oo, (B.7) 

2-+1 



We can however extract the divergence by choosing an irreducible basis for the G-functions 
where the leftmost index is different from 1, except for functions of the form z), e.g. 

G(l, 0; z) = G(l; z) G(0; z) - G{0, 1; z). (B.8) 

In this basis the singularity structure in z = 1 is explicit. Note at this point that the 
pentagon integral does not contain M-functions of the form M(l, . . .), i.e. , I and ^ are 
well-defined. 

In Ref. [28, 39], an algorithm was presented which allows to express 2dHPUs of the 
form H{'w{a); z) in z = 1 in terms of ordinary HPVs in a. First write H{w{a); 1) as the 
integral of the derivative, 

H{w{ay, 1) = H{w{aoy, 1) + £ da' H{w{a')- 1), (B.9) 

where oq is arbitrary, provided that H{w(aQ); 1) exists. The derivative is now carried out 
on the integral representation of -ff, and then we integrate back. Since differentiation 
lowers the weight by one unit and since we know recursively how to reduce all H{w'{a) \ 1) 



'^For the meaning of the indices mi and Zi , the reader should refer to Ref. [36] . They are equivalent to 

the corresponding notation for HPL's, e.g. H{0, 1; z) = H2{z). 

®In some cases the divergence is tamed, e.g. lim^^i Inz ln(l — z) = 0. 
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with w' < w, we obtain now the function H(w{a); 1) with a being the Hmit of the last 
integration. 

Since M-functions are nothing but the values in z = 1 of the G functions, and since 
G-functions are straightforward generalizations of the HPUs, this algorithm immediately 
generalizes to the M-functions, and allows in principle to reduce functions of the form 
M(. . . , Ai, . . .) to HPUs of the form H{. . . ; Ai), the algorithm being a straightforward 
generalization of Eq. (B.9) 

M{w{a)) = M{w{ao)) + [ da' M{w{a')), (B.IO) 

where oq is arbitrary, provided that M{w(ao)) exists. 

C. Generalized hypergeometric functions 
C.l Appell functions 

The Appell functions are defined by the double series, 

Fi{a,b,c,d;x,y) = ^ -— — -, 

m=0 n=0 {(^)m+n m. n. 

oo oo 



„ / , , X (a)m+n ib)m (c)„ x"^ y"" 

F2{a,b,c,d,e;x,y) = > > TTTT^^ TH"' 

f^n (d)m (e)„ ml n! 



m=0 ii=0 
oo oo 
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Prom the series representation the Mellin-Barnes representation is trivially obtained. For 
some special values of the indices, the Appell functions reduce to simpler hypergeometric 
functions, e.g. , 
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The analytic continuation of the Appell function reads 

T{d)T{b-a 
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r(o)r(d-6) 

C.2 Kampe de Feriet functions 

The Kampe de Feriet functions are defined by the series 
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with l<i<p, l<j<Q', l<A;<p' and 1 < i < q', and convergence requires 
p + q < p' + q' + 1 [29] . The order of a Kampe de Feriet function is defined hy uj =p' + q'. If 
p + q = uj + 1, the function is called complete, and all other cases are obtained as a confluent 
limiting case of a complete function. From its definition it is clear that the Kampe de Feriet 
function enjoys the following symmetry property, 
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The Kampe de Feriet functions encompass the Appell functions for particular values of the 
parameters, 
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Fi{a,b,c, d;x,y), 
F2{a,b,c, d, e;x,y), 
Fs{a,b,c, d, e;x,y), 
F4{a,b,c, d;x,y). 
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The Kampe de Feriet function involves only Pochhammer symbols of the form (.)„^+„2, 
(•)ni , (•)n2- We could alternatively define a function involving (.)„j_„2- Let us consider the 
double series 



p',q',r' 



bj Ch 
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Note that this function can appear in the analytic continuation of the Kampe de Feriet 
function. In the following we proof that the generalized hypergeometric series F can always 



-43 - 



be reduced to Kampe de Feriet-type functions. We rewrite the double sum according to 

oo ni 



oo oo 



oo oo 



ni=0 712=0 ni=0 n2=0 ni=0 n2=ni ni=n2=0 

oo ni oo n2 oo 

ni=0 712=0 n2=0 ni=0 ni=n2=0 

The double sums are now reshuffled into double sums up to infinity, e.g. 

oo ni oo oo 

^ ^ ^ ^ ^ , with ni = N + n2, 

ni=0 712=0 N=0 n2=Q 
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and equivalently for J2'^2=o Sni=o- second sum we then rewrite all Pochhammer 

symbols with negative index according to 
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We then arrive at 
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We see that in Eq. (C.ll) only Pochhammer symbols of the form (.)„^+„2, and (.)„2 
appear. Note however that this class of functions encompasses a larger class of functions 
than the one defined by the Kampe de Feriet function. 

The analytic continuation of the Kampe de Feriet function is easily carried out using 
its Mellin-Barnes representation. In particular, we use here the following results, derived 
from Mellin-Barnes integrals. 
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and the F function can reshuffled into Kampe de Feriet-type series using the algorithm 
described in the previous paragraph. 



D. M. functions 

In this section we apply the techniques combining the Mellin-Barnes and series represen- 
tations to study the analytic properties of the M functions that appear in the e-expansion 
of the F4 function, Eq. (5.20). We start by analyzing some general properties of these 
new transcendental functions, and study more specific properties like analytic continuation 
and reduction to known functions in subsequent sections. First, it is easy to see that M 
functions are symmetric, 

M{z,j,k;x,y) =M{3,i,k;y,x). (D.l) 

Second, when one of the arguments is zero, then this function reduces either to or to an 
S-sum at infinity. Indeed, if say y is equal to zero, then only the term n = in the sum 
contributes, and since S^{0) = 0, we get 

• If Jt^ 0, then 

00 / \ 2 

M{iJ,k;x,0) = Y^ i"^] x"-'Sr{m)Sj{0)S,;{m) = 0. (D.2) 

• If J= 0, then 

M{7, 0, k; x,0) = Y^ ( J Sj{m) Sj^{m) = ^ x'" Sj{m) S^{m). (D.3) 

m=0 m=0 

The product of nested harmonic sums can now be reduced using their algebra, and 
so we get 

M{i, 0, k; X, 0) = ^ 5(00; 0, If, X, 1, ... , 1). (D.4) 

i 

where the 5-sums are defined recursively in a similar way as the nested harmonic 
numbers, [36] 

N 

S{N;r;x) = Y,^S{k;^;x'), (D.5) 

k=l 
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and the values of the iS-sums for A'^ = oo are related to Goncharov's multiple poly- 
logarithm. 

We now analyze the analytic continuation of some of these functions under the trans- 
formation y — > 1/y. Due to the symmetry property (D.l) the analytic continuation in x 
follows the same lines. For A; > 1 we have, 

V'(l + n) = -7ij-h5i(n), 

+ n) = k\ (a+i - Sk+i{n)) , 

where 7^ denotes the Euler-Mascheroni constant, and ip^^'^ denote the polygamma functions 

V^^'^l^) = (^^J Inr(z). (D.7) 
We can easily relate the function M. to Mellin-Barnes integrals of the form 

where ^' (2:1,22) denotes any product of polygamma functions with arguments 1 + zi, I + Z2 
or 1 + zi + Z2- For example we can write 
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= M(l,0,0;x,y) -7£;A1(0,0,0;x,y). 

It follows then that studying the analytic continuation properties of certain classes of M.- 
functions is equivalent to study the properties of the McUin-Barnes integral (D.8). If both 
X and y are smaller than 1, then we close both contours to the right, and we take residues, 
and we find the series definition of the A^-functions. 
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Note the appearance of ln(— y), which produces an imaginary part for y > 0. 

In some cases it is possible to express the A^-functions in terms of known functions. 
This becomes possible by relating those functions back to simple cases of the expansion of 
the Appell F4 function, in which cases we can apply the reduction formulae (C.2) given in 
Appendix C. As an example, let us consider the reduction formulas 

FJa,P,a,p- ^ = i^-^ni-y)\ (D.ll) 

V {^-x){l-y) [1 - x)il - y)J 1-xy 
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where x and y are solutions of 



Solving this system we arrive at 



]- (x + y-l + VA(l,x,y)) , j' = ^(^ + y- l + • (D.13) 
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Prom this we can immediately read of the A^-function of weight 0, 



X(0,0,0;x,y) =F4(1, 1,1,1; = ^"^^1, (D.14) 



Using Eq. (5.17), we can easily obtain all the A^-functions of weight 1, e.g. 

d 

M(l, 0, 0; X, y) = — F4(l + e, 1, 1, 1; x, y)|,=o 

= ^ [(1 - 5)1+^(1 - yf+' 2-Fi(l + e, 1 + e, 1; xy)] i^^^ (D.15) 

_ {l-x){l- y) 
l-xy 

Similarly one can obtain 



(ln(l -x) + ln(l -y)-2 ln(l - xy)) . 



A1(0,l,0;x,2/) =iL-^KL_^ (in(l-x)-21n(l-a;y)), 
i xy 

Af(0,0,l;a;,y) = il^^Kl^ (1^(1 - y) - 2 ln(l - xy)) . 
i xy 



(D.16) 



Note however that starting from weight 2 we cannot obtain a reduction in this way in 
all cases, since from Eq. (5.17) we see that by expanding Pochhammer symbols in the 
denominator we only produce <Sii(n), so we cannot obtain S2{n) = 5ii(n) — Zii(n), which 
needs the contribution from a Pochhammer symbol in the numerator. This however cannot 
be achieved by the Appell F4 function, and we have to consider more general functions, 
like the Kampe de Feriet functions. In the particular situation however where the M 
function can be related to an Appell F4 function, we can sum up the series. For functions 
of weight 2 we can sum up all the series except for A^(2, 0, 0; x, y) (and the related function 
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M{0,2,0;x,y) = M{2,0,0;y,x)). 

M ((1, 1), 0, 0; ,) = ii^^i^ ( ^^-^^ + ^^o,\l -y) ^ Wil- .y) 
^ ' 1 — xy) \ 2{x — 1) 4 — Ax I — x 
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Uh n _ (l-ii)(l-^) A2x-3)Li2(xy) , xlog"(l-y) , (x - 2) log^(l - xy) 

jvi^i,u,i,x,y) y 2{x-l) ^ A- Ax ^ x-1 

- 1 log(l - X) log(l - y) + Ml-|)log(l-^y) _ 1 _ J ^ 

7W(0,l,l;x,y) = 7W(l,0, l;y,x), 

(D.18) 

yum n 9- . „^ - _ (.Ty-l)Li2(xy) (4 - 5y) log2(l - x) 



1-xy) \ 2(5;-l)(y-l) 4(y - 1) 

(4-5.x)log^(l-y) (-2xy + x + y)log^(l-xy) 

+ 77^ — 7-\ + F — Tv^ — 2 - a;) log(l - y) 

4(x-l) (x-l)(y-l) 

^ (3y - 2) log(l - x) log(l - xy) ^ (3x - 2) log(l - y) log(l - xy) \ 
y-1 x-1 j' 

E. Analytic continuation of the scalar massless pentagon 
E.l Analytic continuation from Region 11(a) to Region 11(b) 

In this section we explicitly perform the analytic continuation in ti/t2 ^2/^1 (or equiv- 
alently y2 — > l/y2) of the solution (5.13) valid in Region 11(a) to Region n(b), and we 
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explicitly proof the relation (5.15). Let us start with the first term in Eq. (5.13). Using 
Eq. (C.3), we find 



jL (_i)2^ r(l - 2e) r(l + e)2 y^-i F4 (l - 2e, 1 - e, 1 - e, 1 - e; ^) 
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Similarly, the second term becomes, 
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The Kampe de Feriet function have aheady the correct form. Combining the remaining 
Appeh functions, and using the fact that In 2/2 — Inyi = — In and 



-ITT + - (-1)^ r(l - e) r(l + e) = V(l - e) - V'(e), 



we find the desired result 



(E.7) 
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E.2 Symmetry properties in Region I 



In this section we prove that the solution has the correct symmetry properties under 
the exchange of the dimensionless quantities xi and X2, 



(E.9) 



The expression (5.26) apparently violates this relation, due to the explicit appearance of 
the argument Xi/x2 in the hypergeometric functions. In the following we show that if we 
perform the analytic continuation of one of the hypergeometric functions according to the 
prescription xxjx^ ^ x^jxx^ then the resulting expression is manifestly symmetric. 

Since and l^^ are obviously symmetric, we focus here only on the remaining terms. 
Using the Mellin-Barnes representation of the Kampe de Feriet function, we find that 
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3 

We already see in this expression that the terms involving F2 q exchange their roles under 
the transformation ti ^ t2- Let us now show that the same hold true for the remaining 
terms, and let us concentrate on the terms in Eq. (5.29) having a coefficient involving 
(— ii)"*^. The terms involving derivatives of Kampe de Feriet functions become 
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where we used the fact that for any function / 

d 



aj(/W + /M)),,^„ = o. 



(E.12) 



The resulting function is now symmetric to the corresponding one with coefficient (— ^2) ^• 
FinaUy for the terms involving an Appell F4 function we find 



In xi — iiT 



^4(1, 1 — e, 1 — e, 1 + e; —xi, —X2) 



Inxi - - + V'(-e) - V'(e) + (-1)' er(-e) r(e)J ^4(1, 1 - e, 1 - e, 1 + e; -xi, -X2), 

(E.13) 

where the last step follows from Eq. (E.7). This term is symmetric to the corresponding one 
with coefficient {—t2)~^ after the latter has been combined with the first term in Eq. (E.IO), 
and this finishes the proof that Eq. (5.29) is indeed symmetric under the transformation 
{ti ^ t2,Si ^ S2). 

E.3 Analytic continuation from Region 11(a) to Region I 

In this section we show how the solution in Region I, Eq. (5.29) can be obtained by per- 
forming analytic continuation from Region 11(a), Eq. (5.13), according to the prescription 

yi 

Let us start with I^^^""^ and l2^"^ . Using the analytic continuation formulas for the 
Appell F4 function, Eq. (C.3), we find 
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l2\s,Si,S2,ti,t2) 



Performing the analytic continuation for I^^^""^ + 1^^"'* and collecting all the terms we find, 
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Comparing with Eq. (5.29), we sec that we find the correct terms. The analytic continuation 
of /g^^"-* +/g^^"'' follows similar lines. We immediately find the remaining terms in Eq. (5.29). 
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F. Analytic expressions of the integrands 

F.l The analytic expressions of the functions i^^^ and i^^^ 

In this appendix we present the analytic expressions of the functions i^^^ and i^^^ that enter 
the integrals (6.13) and (6.14). 

i(0)(xi,X2)= (F.l) 
-i In^ (l - - f^)) + ln(l - w) In (l - - f^)) + ln(l - v) lnx2+ 

lnt>ln (^1 - t;(^l - f^)) + lni;lnx2 - lna;2 In (^1 - i;(^l - f^)) - i ln2(l - v) - - 
lnt;ln(l -i;) - i In^ X2 - ^. 

iW(xi,X2)= (F.2) 

Li2 ( --\— i^rO - - 1 ) ) - 1^ ^2Li2 ( -pT^rO - 

i ln^(l — t;) In ^ — t;^ + vxi — VX2 + v + X2^ — ^ In^ w In ^ — + vxi — VX2 +v + X2^ — 

\ \V? [l - v[l - j In ^ - i;2 + yxi -VX2+V + X2) - 

^ In^ X2 In ^ — + vxi — VX2 + v + X2^ — In ln(l — i;) In ^ — i;^ + vxi — ■UX2 + v + X2^ + 

ln(l - ■«) In ^1 - ■u^l - 1^^^ In ^ - + vxi — VX2 + v + X2^ + 

ln(l — v) In X2 In ^ — ■y^ + vxi — t;x2 + v + X2^ + 

In v In (|l - w (^1 - ^ In ^ - + t>xi — ■UX2 + v + X2^ + 

In V In X2 In ^ — t)^ + vxi — VX2 + v + X2^ — 

In X2 In ^1 — v ^1 — ^ In ^ — i;^ + t;xi — VX2 + v + X2^ — 

- v'^ + vxi - VX2 +v + X2^ - 21n^(l - i;) In (l - v(^l - f^^^ - | ln^(l - i;)lnx2 + 
ln(l - v) In^ (1 - V (^1 - - Ini; In^ X2 + In^ X2 In (l - ^^(l - f^)) + 5 In^ v lnx2 + 

ilnx2ln2 (l-v(l-^yj -21ni;ln(l-'y)ln(^l-'y(^l- 1^)) - Ini; ln(l - i;) lnx2 + 
ln(l-i;)lnx2ln {l-v(l-^y^ + ^tt^ In (l-v(l-j^^^ -ln'ylnx2ln (l-v(l-j^^^ + 
In^(l-i;)+21n'i;ln2(l-i;) + ln2?;ln(l-t;) + |7r2 ln(l-i;)- ^tt^ Inv+^ln^ X2 + Itt'^ Inx2-C3- 

F.2 The analytic expressions of the functions j'^^^ and j^^^ 

In this appendix we present the analytic expressions of the functions j^^^ and j^^-* that 
enter the integrals (6.49) and (6.50). 

/'^Hyi,y2)= (F.3) 

^In^ [v{y2-l) +1) +ln(l-?;)lnyi - ln(l - ?;) In {y{y2 - l) + l) +lnt;lnyi- 
Invln (y(y2-l) + l) -Inyi In {v{y2-l) + l) + \hi^{l - v) + + hiv\n{l - v) + 

1 1„2 „ I TT^ 

2 in yi + 
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j^'Hyi,y2)= (F.4) 

-Li3(f^) -Ml-^)Li.(^g^) -ln.Li,(fi=i]±i) -ln,,Li.(fi=i^) + 

(f^) In (^^2 - 1) + l) + i ln2(l - In ( - + v{y, + - 1) + l) + 
^In2vln(^-'y2yi+i;(yi+y2-l) + l) + 5 In^ yi \n v'^yi + v{yi + y2 - 1) + + 
\ In^ (^^;(y2 - 1) + l) In ( - v'^yx + v{yx + - 1) + l) + 
In i; ln(l - In ^ - v^y\ + + ^2 - 1) + + 
ln(l - i;) In yi In ( - + + ^3 - 1) + 1 ) - 



ln(l - t;) In (%i{y2 - 1) + l) In ( - v^yx + + y2 - 1) + l) + 
In v\a.y\ In (| - u^yi + + 7/2 - 1) + — 
In t; In ^t;(j/2 — 1) + 1^ In - + f (yi + y2 — 1) + — 
Inyi In {v{y2 - 1) + In ^ - v^yx + i;(yi + y2 - 1) + l) + 

^TT^ In - -y^yi + ^(yi + y2 - 1) + 1) - 2 ln2(l - u) In yi + 2 ln2(l - In (^u(y2 - 1) + l) - 
ln(l - v) It? yi - ln(l - v) In^ (^v{y2 - 1) + 1^ - 2 In w ln(l - v) In yi + 
21ni;ln(l - v) In ('z;(y2 - 1) + l) + 21n(l - v) \nyi In (v{y2 - 1) + 1^ 



ivr^ In (^?;(y2 - 1) + 1 j - ln^{l -v)-2lnv In^il - v) - In^ v ln(l - t;) - Itt^ ln(l - v) + 
1^2 , 1^2 



.7r^lnu + iTT^lnj/i + C: 
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G. Reduction of Li„ in Region I 

In this section we prove the reduction formulae of polylogarithms of the form Li„ ( -7 — - — ^ ) 

to G-functions, where a{x\,X2-,v) = y^J-^T^^-)j^x2 ' ^^^^ polylogarithm as Gon- 

charov's multiple polylogarithm 

— -)= -M (0„_i,a(xi,X2, / ( — ) o Q (a(a;i, X2, u)) . (G.l) 

a{xi,X2,v)J \ ^ Jo \tj 

We can now extract the dependence on the variable v using the reduction algorithm pre- 
sented in Appendix B, 

M(0Vi,a(x,,.2,.)) =-Li„(^^^-^)+£ d.^ (I) ol^n{a{x„X2,v')), 

(G.2) 

and integrating back we find the expression of M (0n-i,a{xi,X2, v)^ in terms of G-functions 
of the form G{. . . ,v). However, we still need to face the problem of how to choose the value 
for vq. The 'best' choice is of course vq = 1, since we know that at this point the G-functions 
collapse to M-functions, i.e. Goncharov's multiple polylogarithm ^ . However, it is easy to 



^We could of course choose another value, e.g. vo = 1/2, but then we need to face the problem how to 
obtain the values G(. . . ; 1/2). 
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see that 

lim Li„ ( — ) = oo. (G.3) 

v-*i \a{xi,X2,v) J 

We therefore choose = 1 — e, and we compute 
M (0n-i,a{xi,X2,vj^ 

= ^ { - "^'^ (a(x„xll-e)) -I i (t) ° Hxux2,v')) | 

"^70 "^"L \T) ° a^^K^l'^2,.'))• 



We will show in the following that in the cases n < 3 this limit is well-defined. 
G.l Reduction of Lii 

For n = 1, we do not need the machinery of the reduction algorithm, but the reduction 
can be performed in a straightforward way: 

LiW ' . l = -ln(l-"^^^-^^^+^^ 



a{xi,X2,v)J \ v{v — l) 

= Inv + ln(l — v) — In {v{l — v) + v{xi — X2) + X2) 

= lnv + ln(l - -u) - In (-(Ai - v){X2 - v)) (G.5) 

= Int; + ln(l -v)-\n (-A1A2) - In ^1 - - In ^1 - ^ 

= G{0;v)+G{l;v)-lnx2-G{Xi;v)-G{X2;v), 

where the last step follows from Eq. (6.19). 

G.2 Reduction of Li2 

For n = 2, let us start by evaluating the second integral in Eq. (G.4). The last integral 
over V involves Lii ( / \ ^ ] , which we know recursively from the previous paragraph. 
We immediately find 



= G(0, 0; v) + G(0, \;v)-G (0, Ai; - G (0, A2; i;) + G(l, 0; v) + G(l, 1; v) (g.6) 
- G (1, Ai; x;) - G (1, \2\v)-G (A3, 0; i;) - G (A3, \\v)^G (A3, \v,v) 
+ G (A3, A2; v) - G(0; v) In X2 - G(l; v) In X2 + G (A3; v) In 0:2. 

We now turn to the limit. The integral is obtained immediately by putting v = 1 — e in 
the previous expression. As we are interested in the singular behavior around f = 1, we 
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switch to the irreducible basis where this singular behavior is explicit (See Appendix B). 
Furthermore, since e is small, we can write 

-Li2 (— ^ = -ln^e-lnxilne+^ln2a;i + ^ + 0(e) 



a{xi,X2,l - e) J 2 2 6 

= G(l, 1; 1 - e) - lnxiG(l; 1 - e) + ^ In^ xi + ^ + 0(e). 




lim <^ 1; 1 - e) - lnxiG(l; 1 - e) - 1; 1 - e) 

e— >0+ I 



I 



- 1 - e) [In (1 - Ai) - In (-Ai) + In (A2 - 1) - In A2 + lnx2] } + {finite) 
(finite), 



(G.8) 

where the last step follows from (1 — 1/Ai)(l — I/A2) = xi/x2- 
G.3 Reduction of Lis 

For n = 3, let us start by evaluating the second integral. The last integral over v involves 
Li2 ( —( — - — Y ) , which we know recursively from the previous paragraph. We immediately 
find 

-iG(0; V) W- XX - \G{\\ v) In^ xi + \G (A3; v) In^ xi + G(0; v) In (^2) Inxi + 

G(l; In (x2) In xi - G (A3; v) In (^2) In xi - G(0; v) \d? X2-G{l-v) In^ X2 + G (A3; v) In^ X2 - 

i7r2G(0; v) - \^''G{1- + ivr^G (A3; v) - G(0; v)G (0, Ai; 1) - G(l; t;)G (0, Ai; 1) + 

G (A3; v) G (0, Ai; 1) - G(0; v)G (0, A2; 1) - G(l; v)G (0, A2; 1) + G (A3; t;) G (0, A2; 1) + 

G(0; v)G (Ai, 1; 1) + G(l; v)G (Ai, 1; 1) - G (A3; v) G (Ai, 1; 1) + G(0; v)G (A2, 1; 1) + 

G(l; v)G (A2, 1; 1) - G (A3; v) G (A2, 1; 1) - G(0; v)G (A3, 0; 1) - G(l; v)G (A3, 0; 1) + 

G (A3; v) G (A3, 0; 1) - G(0; v)G (A3, 1; 1) - G(l; v)G (A3, 1, 1) + G (A3; v) G (A3, 1; 1) + 

G(0; v)G (A3, Ai; 1) + G(l; t;)G (A3, Ai; 1) - G (A3; v) G (A3, Ai; 1) + G(0; v)G (A3, A2; 1) + 

G(l; v)G (A3, A2; 1) - G (A3; t;) G (A3, A2; 1) - G(0, 0, 0; v) - G(0, 0, l-v)+G (0, 0, \i;v) + 

G (0, 0, A2; f ) - G(0, 1,0; z;) - G(0, 1, 1; t;) + G (0, 1, Ai; + G (0, 1, A2; v) + G (0, A3, 0; v) + 

G (0, A3, l;v)-G (0, A3, \i;v)-G (0, A3, A2; v) - G(l, 0, 0; v) - G(l, 0, 1; v) + 

G (1, 0, Ai; 1-) + G (1, 0, A2; t;) - G(l, 1, 0; t;) - G(l, 1, 1; 1-) + G (1, 1, Ai; I-) + G (1, 1, A2; + 

G (1, A3, 0; v) + G (1, A3, 1; - G (1, A3, \i]v) - G (1, A3, A2; v) + G (A3, 0, 0; v) + 

G (A3, 0, l;v)-G (A3, 0, Ai; t;) - G (A3, 0, A2; v) + G (A3, 1, 0; t;) + G (A3, 1, 1; t;) - 

G (A3, 1, Ai; t;) - G (A3, 1, A2; v) - G (A3, A3, 0; t;) - G (A3, A3, 1; t-) + G (A3, A3, \i;v) + 

G (A3, A3, A2; f ) + G(0, 0; In + G(0, 1; 7;) In - G (0, A3; v) In X2 + G(l, 0; v) In X2 + 

G(l, l;f)lnx2 - G (1, A3; lnx2 - G (A3,0;f) lnx2 - G (A3, l;f) lnx2 + G (A3, A3; lnx2. 
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We now turn to the limit. The integral is obtained immediately by putting f = 1 — e in 
the previous expression. As we are interested in the singular behavior around f = 1, we 
switch to the irreducible basis where this singular behavior is explicit. Furthermore, since 
e is small, we can write 

-Li2 — = - - In e + - In xi In e - - In In e - -tt In e + - In xi 

\a[xi,X2,i- — e) J 6 2 2 60 

+ ^vr^ Inxi 
o 

11 1 
= -ln^a;i - -G(l;l - e) In^xi + G(1,1; 1 - e)lna;i + -tt^ Inxi 



-i7r2G(l;l-e)-G(l,l,l;l-e). 



We now find 

= lim < ^ln^a;i - ^G(l;l - e)ln^xi + G(l,l; 1 - e)lna;i + ^TT^lnxi 
e-»o+ 16 2 6 

- Jvr^Gll; 1 - e) - G(l, 1, 1; 1 - e) + G(l, 1, 1; 1 - e) 





(G.IO) 



+ G(l, 1; 1 - e) [-G (Ai; 1 - e) - G (A2; 1 - e) - Inxs] 
+ G(l;l-e) 



G (A3; 1 — e) In 0:2 + - (3 In^ xi — 6 In X2 In xi + 6 In^ X2 + tt^) 



+ (finite) 
= (finite), 

(G.ll) 

where the last step follows from Eq. (6.19), as well as 

G(Ai;l) =ln(^l-i- 

G(A2;1) =ln(^l-j-], (G.12) 
G(A2;1) =ln(^l-^ 

H. Reduction of Li„ in Region II 

In this section wc prove the reduction formulae of polylogarithms of the form Li„ ^ ^^^^ 

to G-functions, where a(xi,X2,v) = ■ Since most of the discussion is exactly the 

same as in Appendix G, we will be brief on this. We can start by writing down an equation 
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similar to Eq. (G.4). We then reduce the function Lii in a straightforward way, and continue 
in a bootstrap to reduce Li2 and Lis. Particular care is again needed when evaluating the 
limit that appears in Eq. (G.4). Let us illustrate the cancellation of the poles. For n = 2, 3 
we can write 

. . . .... . ,^2 

— Lio 



«(2/i,2/2,l - e)y 2 



iln^.-ln(|)ln..iln^(|).^..,) 



lln^(|)-C(l;l-e)ln(|)+^(M;l-e)4 + e.(.). 

(H.l) 



Va(yi>y2,l-e)y 6 {yj 2 \yj 2 \yj 6 



m 



-— ^--TT^ Ine + O € 

D fa 

= i ln3 ( yi) - \g[1- 1 - e) In^ f ^) + 1; 1 - e) In f ^ 

6 \yij 2 Vyi/ Vyi 

1 2i„/2/2\ 1 2. 



+ -TT^ In (^^J - g7r^G(l; 1 - e) - G(l, 1, 1; 1 - e) + 0{e). 

(H.2) 

Inserting these expressions into the limits that appear in Eq. (G.4), one can easily show 
that the limit is finite. 
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